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The analysis of canonical vacuum general relativity by R. Beig and N. O Murchadha (Ann. Phys. 174 
463-498 (1987)) is extended in numerous ways. The weakest possible power-type fall-off conditions for the 
energy-momentum tensor of the matter fields, the metric, the extrinsic curvature, the lapse and the shift are 
determined which, together with the parity conditions, are preserved by the energy- momentum conservation 
law T°-''-i, — and the evolution equations for the geometry. The algebra of the asymptotic Killing vectors, 
defined with respect to a foliation of the spacetime, is shown to be the Lorentz Lie algebra for slow fall-off 
of the metric, but it is the Poincare algebra for 1/r or faster fall-off. 

It is shown that the applicability of the symplectic formalism already requires the 1 /r (or faster) fall-off 
of the metric. The connection between the Poisson algebra of the Beig-0 Murchadha Hamiltonians (and, in 
particular, the constraint algebra) and the asymptotic Killing vectors is clarified. Their Hamiltonian HlK""] 
is shown to be constant in time modulo constraints for those asymptotic Killing vectors K*^ that are defined 
with respect to the foliation by the constant time slices. 

The energy-momentum and angular momentum are defined by the boundary term QlK"'] in H[K°'] even 
in the presence of matter. Although the energy-momentum is well defined even for slightly faster than the 
r~^/^ fall-off, we show that the angular momentum and centre-of-mass are finite only if the metric falls off 
as 1/r or faster. QIK"^] is constant in time for those if 's that are asymptotic Killing vectors with respect 
to the foliation by the constant time slices. If the foliation corresponds to proper time evolution (i.e. its 
lapse tends to 1 at infinity), then Q[K°-] reproduces the ADM energy, the spatial momentum and spatial 
angular momentum, but the centre-of-mass deviates from that of Beig and O Murchadha by the spatial 
momentum times the coordinate time. The spatial angular momentum and the new centre-of-mass form an 
anti-symmetric Lorentz tensor, which transforms in the expected way under Poincare transformations. 



* Dedicated to Jim Nester on the occasion of his 60th birthday. 
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1. Introduction 

The quantum field theoretical investigations of the early sixties showed that, strictly speaking, the observables 
of quantum fields must be associated only with finite but extended spacetime domains, i.e. they are quasi- 
local [1]. Quantities associated with spacetime points are not observables, and the global quantities, e.g. 
the total energy or electric charge, should be considered as the limit of quasi-locally defined quantities. 
Interestingly enough (although by different reasons, but) the situation is very similar in general relativity: 
energy-momentum and angular momentum cannot be associated with the points of the spacetime. Any 
such local expression is necessarily pseudotensorial and/or internal gauge dependent. Thus if we want to 
characterize the gravitational 'field' by observables finer than those associated with the whole (necessarily 
asymptotically flat) spacetime, then these observables must also be defined quasi-locally. 

In the last two decades a lot of efforts was concentrated on the investigations both of the general 
framework in which the quasi-local energy-momentum and angular momentum should be constructed and 
the specific constructions themselves (and their properties). Although there is no consensus at all in the 
relativity community even about general questions e.g. when to consider a specific construction to be 
'reasonable', it is naturally expected that the globally defined observables, e.g. the ADM energy-momentum, 
must be recoverable as an appropriate limit of the corresponding quasi-local quantities. Although the energy- 
momentum, both at the spatial and null infinity, is well understood, the (relativistic) angular momentum 
(especially at the null infinity) needs more investigations. In particular, one should clarify the limit of the 
spatial angular momentum and centre-of-mass of Brown and York [2] and the ones based on Bramson's 
superpotential and the use of the holomorphic/anti-holomorphic spinors [3] at the spatial infinity. 

One of the most elegant introduction of the ADM conserved quantities at the spatial infinity is based 
on the requirement of the differentiability of the Hamiltonian. This approach of Regge and Teitelboim [4] 
was refined later by Beig and O Murchadha [5] , recovering the ADM energy and linear momentum and the 
spatial angular momentum of Regge and Teitelboim, but giving a different expression for the centre-of-mass. 
The recent investigations of Baskaran, Lau and Petrov [6] show that the Brown- York centre-of-mass tends 
to the expression of Beig and O Murchadha. 

The traditional ADM approach of the conserved quantities and the Hamiltonian analysis of general 
relativity is based on the 3-1-1 decomposition of the fields and the geometry. Hence it is not a priori clear 
that the energy and the spatial momentum form a Lorentz vector, or the spatial angular momentum and 
centre-of-mass form an anti-symmetric Lorentz tensor. To ensure the Lorentz covariance of the conserved 
quantities at spatial infinity Nester developed a spacctimc-covariant Hamiltonian formulation of general 
relativity [7]. However, the content and the results of the theory can be spacetime covariant even if its 
form is not. Thus, in particular, we should find a spacetime interpretation of the traditional Hamiltonian 
formulation and the 'conserved' quantities of the theory in terms of some appropriately defined asymptotic 
spacetime Killing vectors. 

It is known that the ADM energy-momentum can be finite and well defined (i.e. independent of the 
background structure) and the energy non-negativity can be proven even if the metric falls off with the radial 
distance r slightly faster than [8-12]. This raises the question of finding the weakest possible fall-off for 
the metric and extrinsic curvature under which the spatial angular momentum and centre-of-mass are still 
finite and well defined. Furthermore, we should be able to treat not only the vacuum theory, but the matter 
fields should also be included. 

The present paper is devoted to the investigations of the global energy-momentum and angular momen- 
tum introduced at the spatial infinity of asymptotically flat spacetimes. We extend the analysis of canonical 
vacuum general relativity by R. Beig and N. O Murchadha [5] in numerous ways: the interpretatableness of 
the results in the spacetime is required, the symplectic structure is considered not to be fundamental and 
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the emphasis is shifted to the field equations, the 1/r and a priori fall-off conditions for the metric 
and the canonical momentum, respectively, are relaxed, the matter fields arc included and the background 
dependence of the angular momentum is investigated. In the literature several mathematically inequivalent 
model for the spatial infinity have been suggested (see e.g. [13-16]). However, the notion of asymptotic 
flatness at the spatial infinity based on a spacelike hypersurface is expected to be the weakest possible in 
the sense that in every reasonable; model of spatial infinity the existence of such a hypc^rsurfacc is expected. 
Thus, in the present paper, we do not use any specific model of infinity, and the notion of asymptotic flatness 
that we use is based on the existence of a certain spacelike hypersurface. Since there is some recent interest 
in higher dimensional (Lorentzian) models (see e.g. [17-18]), and since no extra effort is needed to do the 
analysis in general m = n+1 spacetime dimensions, we assume only that the dimension of the spacetime is 
m > 3. 

In the traditional analysis the lapse and the shift are implicitly assumed to depend only on the spatial 
coordinates, but not on the time coordinate. However, the equations of motion allow their time dependence. 
It turns out that it is precisely this freedom that makes possible to give the spacetime interpretation of the 

Poincarc algebra of the Hamiltonians found by Bcig and O Murchadha. If we excluded the time dependence 
of the lapse and the shift then we would not be able to recover the boost Killing vectors even of the Minkowski 
spacetime. 

The philosophy of our analysis deviates slightly from the traditional one. It is the equations of motion 
that are considered to be fundamental, and the boundary conditions at infinity arc required to be the slowest 
power- type fall-off conditions compatible with the evolution equations. Then the symplcctic structure and 
the Hamiltonian are considered to be only as secondary structures. They are considered to be important only 
from the point of view of finding observables and, in particular, the conserved quantities, but not from the 
point of view of the boundary conditions. This departs from the philosophy of Regge and Teitelboim, where 
the boundary conditions and the Hamiltonian were determined in a single procedure from the regularity and 
differentiability of the Hamiltonian. Thus, having the boundary conditions been specified, the Hamiltonian 
H must be chosen such that the correct field equations be recoverable as the flows corresponding to the 
Hamiltonian vector fields of H. Then the value of this Hamiltonian on the constraint surface will define 
the ADM quantities. Although the Beig-0 Murchadha analysis was carried out for the vacuum Einstein 
theory, the value of their Hamiltonian on the constraint surface can be used to define the energy-momentiim 
and angular momentum even in the presence of matter fields. We clarify in what sense these quantities are 
conserved, and how they depend on the flat background metric. 

First we determine the weakest possible power-type fall-off conditions for the energy-momentum tensor 
of the matter fields, the metric, the extrinsic curvature, the lapse and the shift which, together with the parity 
conditions of Regge and Teitelboim, are preserved by the energy-momentum conservation law T"**;;, = and 
the evolution equations for the geometry. In an n -|- 1 dimensional spacetime they are of order 0(r~*-"^^-'), 
0(r~'^) for some fc > 0, 0{r~^^~^^^), 0{r) and 0(r), respectively. The spacetime vector fields built from 
these allowed lapses and shifts will be called the allowed time axes. Then the asymptotic spacetime Killing 
vectors (with respect to an allowed time axis ^") are defined to be those vector fields for which the Killing 
operator is of order 0{r^^)^ and the space of these asymptotic Killing vectors will be denoted by A^^ . Its 
factor /Q^ by the subspace C Af (whose elements are the asymptotic Killing vectors with 0{r^~'') 
asymptotic behaviour) can be endowed with a Lie algebra structure in a natural way, and this is shown to 
be isomorphic to the Lorentz Lie algebra for slow fall-off k S (0, 1), but for faster fall-off, fc > 1, it is the 
Poincare algebra. Thus the structure of the Lie algebra of the asymptotic symmetries is linked to the fall-off 
rate k of the metric. 

One way of associating conserved quantities to asymptotically flat spacetimes is the use of the sym- 
plectic/Hamiltonian formalism. Since the details of the formalism depend on the type of the fields, we 
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concentrate only on the vacuum theory. It is shown that the apphcabiHty of the symplectic formahsm to 
the vacuum general relativity, in particular the existence of the symplectic 2-form, already implies that 
k > ^{n — 1). This excludes the slow {k < 1) fall-off in spacetime dimensions greater than 3. The constraint 
functions are shown to be finite, functionally differentiable and close to a Poisson algebra precisely for those 
lapses and shifts that correspond to the special allowed time axes with 0(r^~'^) asymptotic behaviour. We 
show that the generators of the gauge transformations, i.e. the functions whose Hamiltonian vector fields 
span the kernel of the pull back to the constraint surface of the symplectic 2-form, are precisely these special 
time axes. A subspace C Af is found such that the Hamiltonian H of Beig and O Murchadha, mapping 
Af into the Poisson algebra of functions, preserves the Lie product of the elements of A'^. It is shown that 
H[K°'] is constant in time with respect to the allowed time axes if K'^ S A^, and it is only constant 
modulo constraint functions for K"- G A^^ . 

The boundary term in the Hamiltonian i/[_ftr"] is used to define the energy- momentum and angular 
momentum even in the presence of matter, independently of any symplectic structure. It is shown that, 
although the energy- momentum is well defined even for k > |(n — 2), the angular momentum and centre- 
of-mass are finite only if fc > i(n — 1). Similar result was obtained independently by Baskaran, Lau and 
Petrov in 3-1-1 dimensions recently [6]. HlK""] reproduces the ADM energy and spatial momentum, the 
spatial angular momentum of Regge and Teitelboim and the centre-of-mass of Beig and O Murchadha for 
K'^ e A^ if the allowed time axis corresponds to a gauge generator. However, the familiar boost Killing 
vectors of the Minkowski spacetime are contained in A^ only if the lapse part of does not tend to zero at 
infinity. For K" E A^ with time axes ^" describing pure time translation at infinity the Hamiltonian H[K°'] 
reproduces the energy-momentum and spatial angular momentum above, but gives an additional term (the 
spatial momentum times the coordinate time) to the centre-of-mass of Beig and O Murchadha. To derive the 
familiar, expected transformation law for the (relativistic) angular momentum this extra term is needed. For 
K"^ e A^ the value of the Hamiltonian H[K°'] is shown to be constant in time with respect to the allowed 
time axis provided the constraint equations are satisfied. We investigate the conditions of the background 
(in-) dependence of the energy-momentum and angular momentum, and we found that although the former 
is well defined even if the diffeomorphisms representing the ambiguity of the background metric tend to rigid 
Euclidean transformations as O(r^), where R < —k and R < (3 — n), the latter is well defined if i? < (1 — fc) 
and R < (2 ~ n) (and in the case of the equality, R = (2 — n) , the generator of the diffcomorphism has odd 
parity). In particular, to have well defined angular momentum in 3-1-1 dimensions the metric must fall off 
at least as 0{r~^), and the allowed diffeomorphisms must tend to rigid Euclidean transformations at least 
as O(r-i). 

In subsection 2.1 the necessary tools are introduced and reviewed, mostly to fix the notations and 
conventions. The new key element here is the n + I decomposition and analysis of the Killing operator. 
To motivate the boundary conditions and the precise definition of the asymptotic spacetime Killing vectors 
we discuss the Minkowski spacetime in subsection 2.2. Then, in subsections 2.3 and 2.4, the boundary 
conditions and the asymptotic Killing vectors are discussed. Section 3. is devoted to the analysis of the 
canonical general relativity, in particular to the constraints, the gauge transformations and the Hamiltonian. 
In section 4. we apply the Beig-0 Murchadha Hamiltonian to define the energy- momentum and angular 
momentum of general asymptotically flat spacetimes, even in the presence of the matter fields, and clarify 
how these quantities depend on the background metric. The appendix is the brief discussion of the boundary 
conditions for the matter fields at the null infinity. 

We use the abstract index formalism, and only the underlined and boldface indices take numerical 
values. The spacetime dimension and the signature will be assumed to be to = n -|- 1 and 1 — n, respectively. 
The Riemann and Ricci tensors and the curvature scalar e.g. of the spacetime connection Va are defined 
by -'"i?%ed^°y^Z<^ := Vr(VzX«) - Vz(VrX«) - V[y,z]X«, := ™ii^„e6 and '"i? := ™i^„65"^ 
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respectively. Thus, Einstein's equations take the form "^Rab — ^"^Rgab + ^9ab = —KTab, and we use the units 
in which c = 1. The terminology and formalism that we follow in the symplcctic/Hamiltonian description of 
general relativity are mostly based on [19-22], and the level of the mathematical rigor we work at in section 
3. corresponds to that of [21,22]. 

2. The TO = n + 1 decomposition 

2.1 The n+1 form of the Einstein equations and the KilHng operator 

Let E/ be a smooth foliation of the m — n + 1 dimensional spacetimc {M,gab) by spacelike hypcrsurfaccs 
and let be its future directed unit normal such that f^Vat is positive. Let Pj^ := 6^ — f^tb, the orthogonal 
projection to Sj. The lapse function N, the extrinsic curvature Xab and the acceleration of the foliation 
arc defined by Nf^Vj := 1, Xab ■= PaP^Vetf = ^LtQab and := f^VJe = ~De{lnN), respectively. 
Here cjab is the induced (negative definite) metric and Da is the corresponding intrinsic Levi-Civita covariant 
derivative. The corresponding Ricmann and Ricci tensors and the curvature scalar will be denoted by R°'bcdi 
Rab and R, respectively. If ^" is any smooth vector field such that ^"Vai = 1 ('evolution vector field' or 
rather 'general time axis'), then has the form = Nf^ + N"' for some vector field N"- = N''P^, the shift 
part of For any vector field and purely spatial tensor field T^^'"^^ let us define the 'time derivative' 

K'::::: ■■= •••nr = np-i^.p-^pI^ ...p/;LtT;-;;; + LNT«-:-^ (2.1.1) 

In particular, for the time derivative of the induced metric we have 

qab = 2Nxab + LN^ab- (2.1.2) 

Let £ai_...am be the spacetime volume m-form. The induced volume n-form and volume element on is 

defined by eai...a„ := i"£aai...a„ and dSt := ■^£ai...a„ = ■\/Md"x, respectively, and hence dv = NdSt'^t. 
(The other convention for the orientation of the submanifolds, which would be slightly more convenient if 
unitary spinors were used [23], is when eai...a„ := £oi...o„ot" = (— )"i'^eoai...o„-) If S C S is a compact n 
dimensional submanifold with smooth boundary <S, v"" its outward directed unit normal in S, then, by the 
negative definiteness of qab, for any vector field X"" tangent to S the Gauss law takes the form DbX^dYi = 
~ §gX''vbdS, where dS := (^.^^-^^y fvf eefa2...a^- 

The conservation of the matter energy-momentum, i.e. T°'';(, = 0, is equivalent to 

fi = N(-D,f + a'^'Xaft - J^fD.N - nx) + ^nM, (2.1.3) 
i*- = Ar(-£)„a^'^ - ^a'^D^N + l^j^D^N - 2i„x"' - Xi') + (2.1-4) 

where := T^Hah, j'' := T^HaP^ and a"^^ := T^^P^P^. The projections of the to dimensional Einstein 
equations, "^Gab + ^gab + i^Tab = 0, are the constraints 

Kc:=rt''('"Gab+A5„f, + Kr,fc) ^-l(^R+{x'-XabX"'''))+^ + 1^1^ = 0, (2.1.5) 

KCa := P^t" ('"Geb + A5e6 + KT,b) = - (d,x\ ' -Dax) + KJa = 0; (2.1.6) 



and the evolution equations 

Xcd = N(-Rcd + 2xceX"d - XXcd) + ^NXcd - DcDdN+ 

2 / I 1 X (2-1-7) 

+ (^^^^^^'^ + ^^[-^^^ + l^f^^"' + (^^^^^-^j- 

To check whether the evolution equations (2.1.2) and (2.2.7) preserve the constraints, take the time derivative 
of c and Ca and use (2.1.2)-(2.1.7). We get 

c = -2c''DaN - NDac" + Lnc - 2Nxc, (2.1.8) 

Ca = 2cDaN + NDaC + L^Ca - NxCa- (2.1.9) 

Therefore, if the constraints (2.1.5), (2.1.6) are satisfied at t = 0, then any of their derivatives also vanish, 
and hence the constraints are preserved by the evolution equations (2.1.2), (2.1.7). Since in the present 
paper we are interested in asymptotically flat spacetimes, the cosmological constant A will be assumed to 
be zero. 

If if" =: Mt"" + M'^ is any smooth vector field, where = P^M^, then the n + 1 decomposition of 
the so-called Killing operator, acting on the spacetime 1-form field Ka, is 



Nf't^Vi^aKb) = M + M"DaN - N^DaM, (2.1.10) 
2NPyV^i,Kc) = Ma + [NDaM - MDaN^ - 2NxabM'' - LnM„, (2.1.11) 
P^P^V(,Kd) = D(,M,) + Mxab. (2.1.12) 

Clearly, while the space-space projection of the Killing operator is well defined even on a single hypersurface 
S, the first two projections are well defined only if a foliation St of M , i.e. a lapse function N on E, is 
fixed, and, in addition, their right hand side needs a choice for the shift vector A''" too. Obviously, in a 
generic spacetime the Killing equation V(a-ftrb) = has only the trivial solution. However, f^t^V (a^h) = 
and P^f'^V(biV'c) = can always be solved, i.e. the initial value problem for the system 



M = -M''DaN + N''DaM, (2.1.13) 
Ma = - (nD^M - MD^n) + 2NxabM'' + LnM„, (2.1.14) 

is unconstrained, and the initial value problem for M and M" always has a solution. If = Mf^ + is 
another spacetime vector field, then the n + 1 decomposition of the Lie-bracket of if and K"^ is 



K,K 



+r f LmM - j + M, M + [md°-m - md^m) . 

Thus in this decomposition only the time-time and the time-space parts of the Killing operator appear, but 
not the space-space parts. Therefore, if both K°- and K°- satisfied (2.1.13-14), then the projections of the 
Killing operator on the right hand side of (2.1.15) would be vanishing. 



2.2 Boundary conditions I.: Matter fields in Minkowski spacetime 
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Let (M,gab) be the Minkowski spacetime and K. the Lie algebra of its Kilhng vectors. As is well known, it 
contains an m dimensional commutative ideal Z, consisting of the constant vector fields on M and inheriting 
a natural Lorentzian vector space structure too, and K/I w so(l,n). Fixing a Cartesian coordinate system 
X- = (r, X' ), a = 0, 1, n and i = 1, 2, n, (i.e. adapting the coordinates X- to an orthonormal basis 
of the space of the constant vector fields), the translation and boost-rotation Killing 1-forms are well known 
to take the form iff = V eX^^ and Kf- = X^Ve^^ -X^VeX^^, respectively. Therefore, any Kilhng 1-form 
can be written in the form = RabKf- +TaKi = (X' Ve^J - X^ VgX' ) + 2Bi (X' VeT - rVeX' ) + 
Ti VgX' + TVeT, where Rat = —Rba and are constant, and we used the notations Bj := Rjo and 
T := To. 

Let be a r = const hyperplane. R^ := (^ijX'XJ, and for some i? > let Bn be the solid closed 
ball of E,- with radius R and Sr :— dBji its boundary. If is its outward directed unit normal, then 
1 = v°-DaR = —^Qij v'^DaX^ , where qij are the components of the induced flat metric Qab on E^ in the 
coordinate system {X' }, and hence = If / = /(r, i?, ^) is any function, then let us define its even 
and odd parity parts, respectively, by ^/(r, R, ^) := |(/(t, R, ^) ± /(r, R, — ^))- Let Tq be the future 
pointing unit timelike normal to E^, and define the quasi-local energy, spatial momentum, spatial angTilar 
momentum and centre-of-mass of the matter fields in the n-ball Br, respectively, by taking the flux integral 
of the conserved current KaT"-^: 



Er := j ifOr^VfcdE = ^ (^j iidS^R'"~^dR' = (^j + ndS^R'^-'^dR' , (2.2.1.a) 

Pr -^J <r°V6dS = ^ {^jfDa.X^dS)R'"'~^dR' = j^ {^j+fdS^R'"'~^dR', (2.2.1.6) 

:= j Kl^T^^ndY. = £ (^j 2fv^^ DaX^^dS^R'^'dR' = -2 (^j - j^W^ dS^R'^'dR' {2.2.1. c) 

jio := j Kl^T'^^ndE = {^j /iv' dS^ R'"- dR' - tP' = {^j - nv' dS^R'^^dR' - tP\ {2.2.1.d) 

where d<S is the area element on the unit sphere S. (Strictly speaking, the traditional [non-conserved non- 
relativistic] centre-of-mass is J'° -I-tP' .) Since KaT"^^ is divergence-free for Killing vectors, these quasi-local 
quantities are, in fact, associated with the (n — l)-surface Sr and depend only on Ka'. if E is any compact 
spacelike hypersurface whose boundary 9E coincides with Sr, then the flux integral of T'^'^K^ on E will be 
that on Br. The necessary and sufficient condition of the existence of the R ^ oo limit of these integrals, 
respectively, is^ 



R (f ndSR = R'' (f +Atd5 = o(i?-°), (2.2.2.a) 

R <j> fDaX' dSR = R^ <j> +f dS = o(i?-°), (2.2.2.6) 
Jsr Js 

2R^ I fv^' DaX^'^dSR = -2iZ("+i) / -jlVld^ = o(i?-°), (2.2.2.c) 
Jsr Js 

2R^ i /it;'d5ij = 2i?("+i) / -/xuM5 = o(i?-°). (2.2.2.d) 
Jsr Js 

^ A function /(r) will be called of order o{r^'') if limr^oo(/('')^'^) = 0) and will be called of order 0(r^'^) 
if limr^oo (/('')''*') exists. In particular, o(r+°) will denote logarithmic divergence and o(r~°) logarithmic 
fall-off, while 0(1) means that /(r) tends to a constant at infinity. 



7 



These global integral conditions can be ensured by the explicit fall-off and parity conditions 



^(t.R.^) = V""(r,^) +o(R-"'). (2.2.3.0) 

where f{T, R, is an arbitrary function with even parity, n^"^^ and contribute only to the energy 

and the spatial momentum but not to the angular momentum and centre-of-mass, hence we may call them 
the ADM mass aspect of T"*". Repeating this analysis on boosted hyperplanes S^, := {r' := TCOsh/3 + 
X^tti sinh/3 = const}, where /? G R and q^^aiaj = —1, we obtain that ^ f{T,R, ^) = and, in addition 
to (2.2.3a), 



i' (t, R^^) = ^ ^f^'^Kr, ^) + o(ii— ), (2.2.3.6) 
^' ^) = ^ ) + (2.2.3.C) 

Note that although (2.2.2) could be ensured only by fall-off conditions that are strictly faster than those 
in (2.2.3), but with these conditions we would exclude e.g. the electromagnetic field from our investiga- 
tions, where the typical fall-off of the energy-momentum tensor is R~"^. Obviously, the fall-off and parity 
conditions (2.2.3) are only sufficient, and the global integral conditions (2.2.2) can be satisfied without the 
parity conditions too. The advantage of the fall-off and parity conditions is that they can be given explicitly. 
However, the price that we had to pay for this is that we excluded those field configurations from our inves- 
tigations that satisfy the global integral conditions (2.2.2) but not the explicit fall-off and parity conditions 
(2.2.3). 

One can carry out a similar analysis of the fall-off and global integral conditions that can ensure the 
finiteness of the global energy-momentum and (relativistic) angular momentum at the future null infinity. 
Since, however, in the present paper primarily we are interested in the kinematical quantities defined at the 
spatial infinity, the null infinity case will be discussed only in the Appendix. 

In the rest of this subsection we discuss the conditions under which the ^ni{m + 1) spacetime Killing 
vectors can be recovered, at least asymptotically, from quantities defined on a general asymptotically flat 
spacelike hypersurfacc in the Minkowski spacetime. Since, however, many parts of the following discussion 
are well known from various sources, we sketch only the main points of the argumentation. 

The global Cartesian coordinate system X- = (r, X' ) defines a foliation of the Minkowski spacetime 
by the hyperplanes and gives the 'time axis' (^)", i-e. the corresponding lapse is one and the shift 
is zero. Thus our aim is to determine those conditions under which a coordinate system {t,x'-), based on 
a more general spacclike hypersurfacc S, 'approaches asymptotically' the Cartesian coordinate system 'at 
infinity'. However, to do 'a coordinate system approaches a Cartesian coordinate system asymptotically at 
infinity' to be meaningful we need to use some model of the spacelike infinity of the Minkowski spacetime. 
We choose the classical conformal boundary of Penrose. Thus let (M, gab) be the conformally compactified 
Minkowski spacetime (e.g. the closure of the conformally embedded Minkowski spacetime in the Einstein 
universe, see e.g. [24]), Vl the conformal factor on M such that gab\M = ^"^Im gab and i° G M the point 
on the conformal boundary of (M, gab) representing the spatial infinity. Then the family of hyperplanes 
uniquely determines a family S^- of smooth Cauchy surfaces in M such that i'^ € Sr for all r G R, and 
the future directed ^a(,-unit normal fa of all the surfaces S,- coincide at i°. The compactification of J^r to 
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Tt-r can also be done by the standard inversion transformation: on some open neighbourhood Vr of i*^ in 
tr let X' := R-'^X\ Then the coordinates X' can be extended from K - {i°} to Vr by X' (i") = 0, 
and E? := is a ^afe-orthonormal spatial basis at i" orthogonal to fa- This basis turns out to be 

independent of r. 

Let S be any smooth spacelike hypersurface in M through i°. We say that S := S — {i°} is approaching 

the leaves of the foliation Tit at infinity if the future directed ^afc-unit normal ta of S coincides with fa at 
i". We assume that E is such a hypersurface, otherwise it is called asymptotically boosted with respect to 
the leaves S^. (Note that we have to assume the smoothness of E even at i°, because otherwise its normal 
would not be well defined.) Let ta be the future directed gab-unit normal to S, and let Qat and Qab be the 
induced metrics and Xab and Xab the extrinsic curvatures of S and S, respectively. Then, by the construction, 
{T,, i'^ , Qab; Xab) IS an asymptote for (T,, cjabj Xab) in the sense of [13] (see also [25]), and, in addition, the 
normal directional derivative of the conformal factor, ft t°Vaf^|s, has a extension to E such that 
n(iO) = 0, (L'afi)(i") = 0, (^a-D6ii)(i°) = and (i^^ AZ)eO)(i°) = 0, and Xabk = ^Xab + ^Qab holds. Let 
{x^ } be a local coordinate system on some open neighbourhood U of i° in S such that (i°) = 0, and in 
these coordinates (^|r)"|io = E? and (i") = hold (i.e., in particular, qij (i°) = —Sij holds, and {x^ } is a 
normal coordinate system with origin i" G E). Let := Sij ,t' ,fJ , and define the new coordinates .t' := f~'^x^ 
and radial coordinate distance := ^ij x' = on [/ :— U — {i°}. (In general the coordinates x' are not 
the restrictions to S of the Cartesian coordinates .) The properties of r2||> and tl stated above imply that 
in these coordinates f2 = r~^(l + 0{r~^)) and O = 0(r~^^+^^) for some positive k and h. (In fact, even for 
general asymptotically flat spacetimes when E is not a smooth hypersurface at i*^, one has k,h> 1.) Then 
for the metrics and the extrinsic curvatures we have 

^ijdx'dar" =n-^qijdx' dx^ = + ^q[f + 0{r-^) + o(r-'=))dx' da^ , 

Xijdx'daH =(n-^xij -n-3f2gij)dx'dxj = [■;^x[]'^^^ + o{r-^^+''^)^dx' dx^ 

for some qi'f^ and Xi j^''^ depending only on ^ . Actually, for the Minkowski spacetime, both k and h must be 

greater than or equal to 2. Therefore, {.x' } is an 'asymptotically Cartesian' coordinate system with respect 
to q„t,. and, on U, it defines a (negative definite) flat metric o^ab with respect to which {a;'} is Cartesian. 

To complete {x' } to a spacetime coordinate system {t, x' } (at least on an open neighbourhood of 
[/ C E in M), we need a whole family E^ of such hypersurfaces providing a foliation of this neighbourhood. 
However, to ensure that this spacetime coordinate system (and not only {x' } on the single hypersurface 
Eq = E) approaches the Cartesian one, all the leaves Ei of the foliation must approach the leaves E,- at 
infinity, i.e. the future directed gab-unit normal ta{t) of Ej must coincide with fa at i° for all t G R. 
(Otherwise the foliation E^ would be asymptotically accelerating at i" rather than being inertial.) Then 
there is a natural extension of the spatial coordinates from S to all the leaves Et via the construction above: 
on Et let {a;' } be the inversion of the normal coordinates a;' on St based on the basis {E?} at i° and with 
the origin at i*^. Denoting the tangent of the curves 7(t) := {t, x' ), x' = const., by and defining the lapses 
and shifts by = Nt" + N'' = Nt + N", we get that N = Q-^iV = 0{r^'P) and N" = N'' = 0{r-P), 
where p > 2, because the leaves Et of the foliation are tangent to each other at i°. Therefore, we can write 
Ar = Af(0)(^) + O(r-i). 

Since the leaves both of the foliations St and S,- are tangent to each other at i" and both X' and 
x' are normal coordinates based on the same basis E^ at i°, X' = x' + 0(P), where p > 2. But since 
i?2 = r2(l+0(r-i)) we have X' = x'+0{r^-P), implying that (^fr)" = (afr)"+0(ri-P). Furthermore, the 
angle between the normals fa and ia of S^ and Sj is of order P"^ if the leaves of the foliations are tangent to 
each other in the (p — l)th order. Thus the (hyperbolic) cosine of this angle is taT" = i"'f^gab = 1 + 0(P~^), 
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implying that (^)° = = Nf + N"^ = Ni-^^)"" + 0{r^-'^P) + Oir-P). Therefore, the coordinate basis vectors 
(g^r)" tend to (gyr)° asymptotically, but (^)" tends to (757)" only if p = 2, and the lapse has the form 
A'' = l+0(r~^) (whenever r = t+c+0{r~^), where c is a constant). Expanding the 1-form basis (Vgr, VgX^ ) 
in terms of {te,DgX^), where Dg is the induced derivative operator in the leaves S^, we can rewrite the 
general Killing vector of the Minkowski spacetime given in the first paragraph of this subsection. We obtain 
Ke = Rij {x' Dexi - xi Dex' ) + 2Bi (.'/;' te - tDeX' ) + .Si Dex' + ste, where s{t, x^ ) = s'^'^^t, ^) + O(r-i) 
and Si (t, x^ ) = sp''(t, ^) + 0{r^^). Its structure is similar to that of Ke in the Cartesian coordinates, but 
instead of the constant components of the translations, s and Si are functions of t, ^ and higher powers 
of r~^. Thus they analogous to the supertranslations in the BMS group of the null infinity, and hence it is 
natural to call st" and s' (gfr)" temporal and spatial supertranslations, respectively. 

To summarize: although the global energy-momentum of the matter fields in Minkowski spacetime can 
be ensured to be finite by the i?,^™ fall-off conditions, both at the spatial and the null infinity, to have finite 
angular momentum and centre-of-mass additional global integral conditions must also be imposed on the 
mass aspect of T"*". At the spatial infinity these global integral conditions can be ensured by explicit parity 
conditions. To be able to recover the familiar Killing vectors in their usual form on a general asymptotically 
flat spacelike hypersurface S, at least asymptotically, the lapse N, defining the time coordinate t, must tend 
to 1 as r ^ 00. 



2.3 Boundary conditions II.: Asymptotically flat spacetimes 

Suppose that T, is asymptotically Euclidean in the sense that for some compact subset if C S the complement 
T, — K is diffeomorphic to a finite disjoint union of manifolds S(j), each of which is diffeomorphic to R" — B, 
where S is a solid ball in R". The pieces E(i) are called the asymptotic ends of S. Since the next analysis 
can be repeated on each E(j), for the sake of simplicity we assume that there is only one such end. Suppose 
that there is a (negative definite) metric o9ab on S such that it is flat on the asymptotic end S — . 
Let {a;'} be a coordinate system on — K which is Cartesian with respect to oQab, ^'^ ■= Sijx'-x^ , the 
radial distance function with respect to 0^065 and let o-Dg be the Levi-Civita covariant derivative operator 
corresponding to oQab- Then the quotients ^ can be interpreted as coordinates both on the unit sphere 
S w S"~^ and the sphere Sr of large coordinate radius r in T, — K, and the components qw' of the outward 
directed o9ah-unit normal qI'" to Sr in the coordinate system {x^} too. By a ball of radius r in S we mean 
Br:={peJ:-K\r{p)<r}LlK. 

Let us consider the first, intuitively obvious condition of asymptotic flatness on the components of the 
metric and extrinsic curvature in the coordinate system {a;' }: for some positive k and / 



Qii {^'^ ) = 0* j + ^'^^ ( — ) + o(r-'), (2.3.1a) 

Xij(^'^) = 7Xij('H^)+o(r--'). (2.3.16) 

Therefore, the coefficients gij Xij can be interpreted as functions defined only on the unit sphere S. 

Following [4] and [5], in addition to the fall-off conditions we impose the following global parity conditions 
on the leading terms of the metric and extrinsic curvature: 

,ij W(-^) = W(^), XiJ^')(-^) = -XiS^H^y, (2.3.2a,6) 
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i.e. q-ij ^'^^ is of even, while Xij is of odd parity. Furthermore, we assume that the 'rests' iriab '■= Qab 
oQab — r~''qab^''^ and kab '■= Xab — r~''Xab^^^ satisiy the additional conditions 



oDcmab = o{r ^ oDdoDcTUab = o{r ^ 2), oDeoDdoDcniab = o{r ^),... (2.3.3a) 

oDekab = o{r-^-^), oDdoDckab = o{r-^-^),... (2.3.36) 

which imply that oDe^.-oDe^qab = 0(r"(''+*) ), s = 1,2,..., and, similarly, o-De, ■•■o-DeiXab = 0(r"('+')). 
These properties make the calculations easier. The parity of these derivatives is (— )* and (— )*'''^, respectively. 
The properties rUab = o{r-^), oDcniab = o{r~'^~^), oDdoDcniah = o{r-''-^), oA;,, ■••o-Dei"iab = o(r"(''+^)) 
of the rest niab will be denoted by rUab = o^{r~^), and, similarly, kab = o*(r~'). Although in the actual 
calculations we will use these (essentially technical) assumptions for some finite value of s, for the sake of 
simplicity we assume that rUab = o°°(r~'^) and kab = o°°{r~^). We may call the asymptotic end {'S, Qab, Xab) 
to be (fc, /)-asymptotically flat if for some background metric o9o6 the conditions (2.3.1)-(2.3.3) are satis- 
fied. However, we will see in subsection 4.2 that this notion depends on the background metric too, thus 
(S, Qab) Xab) wiU be Called {k, I) -asymptotically flat with respect to the background metric oQab if (2.3.1)-(2.3.3) 
are satisfied. 

Now let us suppose that /x, j" and a"'^ satisfy the fall-off and parity conditions (2.2.3) with respect to 
the Cartesian coordinates on S with the additional (technical) requirement that the 'rests' are of order 
o°°{r~'^). Next ask what conditions should we impose to ensure the existence of the limit 

Q„[M,M«] := lim / f/xM + i»M„)dS (2.3.4) 
and, considering S to be a leave of the foliation defined by a lapse function N , its time derivative 

Qm[M,M»]= lim / |/i(M- Ar«£)„M + M«£)„iV) + 

+f {Ma - MDaN + NDaM - LnM„ - 2xabM''iV) + ^2.3.5) 

+a''b^{MXab + D(^aMb)) + 

+Da ((AiM + fMb)N'' - {fM + a''^Mb)N^ }( 



f-dS 

with respect to the time axis Nf" + TV". To obtain (2.3.5) we used (2.1.2)-(2.1.4). Note that if 

K'' = Mt" + M", then Q„^[M, M"] is just the n + 1 form of QmiK"] := KoT^HbdY., and Qm[M, M"^] is 
the n + 1 form of f^{T''\VaKb)N + Da{{P^tc - P^tb)i''T'='^Kd))dT.. Thus let us suppose that M and Mi , 
where the latter is defined by Ma =: Mi Do a;' , have the asymptotic form^ 



Mit,x^)=r^M^^Ht,—) + o^{r^), (2.3.6a) 



Mi{t,x^) = r^ Ml^\t,^) + o°°{r^) (2.3.66) 
for some A, B. Substituting these into (2.3.4) we obtain that Qm[M, M"] exists precisely if 



^ To be consistent with our previous notations, we would have to write M*^ instead of M'^^\ However, 
we resolve this apparent inconsistency with the convention that the power of r is always a capital, while the 
power of 1/r is always a lower case letter. In particular, the leading term in the expansion of /(r, ^) can 
be written as r^/^^H^) or r-"/(«)(^). 
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A<1, B<1; (2.3.7) 

and if the equality holds in these inequalities, then the leading terms, M^^'^{t, ^) and M^^^\t, ^), must be 
odd parity Junctions of respectively. The existence of Qm[M, M"] restricts TV and TV". Before discussing 
these restrictions, clarify first the conditions for N and iV" by means of which the conservation equations 

(2.1.3) and (2.1.4), and the consequence (2.1.2) of the definition of the 'time derivative' preserve the fall-off 
and parity conditions (2.2.3) for the matter fields and (2.3.1a-3a) for the metric, respectively. Writing TV and 
A^i in the form of M and Mi given by (2.3.6) with some powers C and D, substituting them into (2.1.3), 

(2.1.4) and (2.1.2) and requiring their right hand side to have 0(r^™) order and even parity, 0(r~™) order 
and even parity and 0(r^'^) order and even parity, respectively, we obtain that 

N(t,x^) = N^^\t,—) + o-=^(r^'), (2.3.8a) 
r 

N, it, x"^ ) = 2x'^puiit) + n it) + vf\t, ^) + o°°(r^), (2.3.86) 

where n [t) and pij [t) = —pj j (t) arc independent of the coordinates {x^ } but may be arbitrary functions 
of the coordinate time t, the powers C and F satisfy 



C<min{l,/-fc}, F<{l-k), (2.3.9) 

and if the equality holds then the leading terms N^'-^\t, ^) and i^i^\t, ^) must be odd parity functions 
of respectively. (The first two terms on the right of (2.3.8b) together is just the kernel of the Killing 
operator o£'(i-/Vj-) = in (2.1.2).) Note, first, that by (2.3.9) F < 1, because we assumed that fc > 0, and, 
second, that there is no reason to keep the term n (t) in (2.3.8b) if > 0. Substituting (2.3.6) and (2.3.8) 
into (2.3.5) and taking into account the conditions (2.3.7) and (2.3.9), one can check that the integral exists 
(and, in particular, the total divergence gives zero). If K"^ = Mt"" + is a spacetime Killing vector then, 
by (2.1.10-12), QmlK""] is zero, as it must be since then T°-''Ki, is divergence-free. 

The notion of the {k, ?)-asymptotic flatness is referring only to the asymptotic end of a single spacelike 
hypersurface S. To ensure that the spacetime itself, i.e. the evolution of S, is also asymptotically flat, we 
must ensure the compatibility of the boundary conditions with the evolution equations for the geometry. 
Thus let us consider the evolution equation (2.1.7) for the extrinsic curvature and ask under what additional 
conditions for N and Na do these equations preserve the fall-off and parity conditions (2.3.1b-3b) for the 
extrinsic curvature. An analysis similar to that we did above yields that (2.1.7) preserves these asymptotic 
conditions precisely when 

N(t, x^ ) = 2x^Pk (t) + T(t) + r^u^^^ it, —) + o°°(r^), (2.3.10a) 

r 

Ni {t, x'') = 2a;''pki (t) + n [t] + r^Vi ^^\t, ^) + o°°(r^), (2.3.106) 
where T{t) and f3i (t) are independent of {x^ } but may depend on t, the powers E and F satisfy 



E<{2-1), F<{l-k), (2.3.11) 

and if the equality holds in (2.3.11) then and fi^^^(^) are of odd parity, respectively. Then by 

(2.3.9) E <C <{l-k), implying that E < {1 - ^k) < 1 and k + E < I < 2- E. In addition, if T{t) ^ then 
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I < k + 2, and, by (2.3.9), k < I even for £; < 0. If A (t) ^ then / = fc + 1, and hence E <(\-k) <\. In 
the presence of matter T{t) also implies I < w and Pi (t) ^ also implies / < w; — 1, where w is the actual 
order of the leading term of the spatial stress (and of the energy density and momentum density), which, as 
we saw, must satisfy w > m. (Interestingly enough, it is just the fall-off conditions k>0, l = k + l that 
ensure the existence of the spinor Chern-Simons functional Y on spacelike hypersurfaces in 3+1 dimensional 
spacetimes [26], by means of which the vacuum Einstein equations can be recovered as the necessary and 
sufficient condition of the invariance of Y with respect to infinitesimal spacetime conformal rescalings [23].) 
For n = 3 and the a priori powers I = k + 1 = 2 the expression (2.3.10) is almost the condition of Beig 
and 6 Murchadha obtained from the investigation of the vacuum evolution equations. The only (and, as we 
will see, important) difference is that the evolution equations, even in the presence of matter, allow the time 
dependence of the coefficients in N and A^". 

2.4 Allowed time axes and the asymptotic Killing vectors 

In the previous subsections the time axis with respect to which the time evolution was defined and the 

generator K'^ of the physical quantity in QmlK""] were treated separately: While the role of was to 
provide a differential topological backgromid, e.g. a foliation of the spacetime and a shift vector to carry 
out the analysis (general time axis) , the nature of K"^ told us whether the corresponding Qm [K'^] should be 
interpreted e.g. as the energy or a component of the spatial angular momentum of the matter fields. Indeed, 
the lapse and shift parts of K" were defined with respect to the foliation that defined. 

However, the structure (2.3.10-11) of TV and N°- is compatible with (2.3.6-7), or, in other words, the 
time axes ^" can be considered as special generators X". (Note that the vector fields ^" and K'^ obtained 
here might be tangent to S or even vanishing, and their lapse part might be positive on some subset and 
negative on other subsets of S. In spite of this fact we call an 'allowed time axis'.) Furthermore, these 
two roles are mixed in the Hamiltonian formulation of the dynamics of the mattcr-|-gravity systems: M and 
M" in the total Hamiltonian (in particular in its matter part Qm[M, M"]) play the role of the n+1 form of 
the generator K"^, and, at the same time, the Hamiltonian generates the time evolution of the states with 
respect to the spacetime vector field = Mt" + M". Thus we assume that all the generators M, Mi and 
M, Mi also have the structure (2.3.10-11) of the allowed time axes and we write 

M = 2a;'' {t) + T{t) + r^/|(^) + o°° (r^) , Mi = 2a;'' iik i (*) + Ti (i) + r^fn + o°° (r^) , 

M = 2a;'' Bk (i) + f{t) + r'^p,^'^^ + o°°{r'^), M; = 2a;''^ki (i) + Ti (i) + r"Pi + o°°{r"), (2.4.1) 

N = 2a;'' /3k (t) + T{t) + r'^v^'^^ + o°°(r-^), = 2a;''pki (*) + n (t) + r^i^i + o°°(r^), 

where E, F, G, H,K,L < (1 — k) and, in the case of equality here, the corresponding coefficient has odd 
parity. The space of the pairs (M, M") on S given by (2.4.1) will be denoted by A. Considering M to be the 

lapse of a (maybe degenerate) foliation of a neighbourhood of S in the spacetime, A can also be interpreted 
as the space of the spacetime vector fields := Mt"" + M"-, where is the future pointing unit normal to 
the leaves of the foliation. 

According to the double role of the components given by (2.4.1), we form another space of spacetime 
vector fields. Namely, if a lapse N is given on E, then let An be the space of those spacetime vector fields 
X" := Mt" -I- that are defined with respect to the (maybe degenerate) foliation determined by S and 
A^. (If N has zeros, then not every element (M, M") of A determines a spacetime vector field in An- if the 
leaves and S^/ of the foliation intersect each other at some point p, then only those elements (M, M") of A 
can determine spacetime vector fields in An, for which M{t,p)tp + M°-{t,p) = M{t',p)t'p + M°-{t' ,p), where 
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tp and tp are the future pointing unit normal of and Sf, respectively, at p.) The structure of the leading 
two terms e.g. of M and Mi resembles to that of the n + 1 decomposition of the familiar spacctime Killing 
vectors of the Minkowski spacetime with respect to a spacelike hypersurface. However, although the first 
terms are linear in, and the second terms are independent of the spatial coordinates, the third terms (which 
would be analogous to the supertranslations of subsection 2.2) may depend on the spatial coordinates and 
may even be diverging. Moreover, although by (2.1.12) P^P^V (cKd) = 0{r~''), i.e. tends to zero at infinity, 
in general neither P^t^V(ftiirc) nor t'^t'^V (aKb) tend to zero. The space An does not form a Lie algebra. In 
fact, if := Mf" + and := Mt" + M'' are any two elements of An, then the first two terms of 
their Lie bracket in (2.1.15) are given by (2.1.10) and (2.1.11), respectively. Substituting (2.4.1) into these 
formulae we find that their asymptotic structure is dominated by N~^x"^ , which deviates from that of 
the allowed time axes. Moreover, while the time dependence of the Killing vectors in a coordinate system 
adapted to the translations is very specific, the time dependence of M and Mi is not specified at all. In this 
subsection we use this freedom to specify a class of spacetime vector fields, whose elements can naturally 
be interpreted as asymptotic Killing vectors. The components of the spacetime vector fields that are only 
allowed time axes will be denoted by Greek letters, as those of and TV" in (2.4.1). 

We noted in subsection 2.1 that the parts P^t^V(tifc) and f^i^V (^aKb) of the Killing operator acting 
on some K"^ G ^jv can be required to be zero. Thus let us define the allowed time axis K"^ = Mt"' + M°' 
to be a strong asymptotic Killing vector with respect to the foliation characterized by the lapse N if its 
components M and M" satisfy (2.1.13) and (2.1.14). However, to ensure that the components of these 
asymptotic Killing vectors be well defined, a choice for the shift vector N"" must also be made. The set of 
these asymptotic Killing vector fields will be denoted by y^^, where i^"' := Nf^ + TV", which will be assumed 
to be an allowed time axis. Obviously, this notion of asymptotic Killing vectors depends sharply on the lapse 
N and the hypersurface S (on which the foliation is based): if a neighbourhood of S in the spacetime is 
foliated by another lapse N from (maybe) another hypersurface E with normal t", then the corresponding 
parts t^i'^V QjKc) and Pji^i'^W (hK^) of the Killing operator will not be zero. (Perhaps, the notation A*^ is not 
very fortunate, and it would have to be denoted by A^^ n)-n<' *° stress that this notion of asymptotic Killing 
vectors depends both on S and N, and the asymptotic Killing vectors themselves are parameterized by 
using the shift N"' too.) Therefore, this notion of the asymptotic Killing vectors appears to be unnecessarily 
strong, and it could be enough to require that the parts Pj^fV^hKc) and f^t^Wf^aKb) of the Killing operator 
be at most of order 0{r~'') asymptotically, i.e. explicitly 

f't^V^aKb) = j^{M + M'^DaN - N-DaM) = r^K(^) (t, ^) + (r^) , (2.4.2) 

2P^t'\/^bKe) = ^ (Ma + {NDaM - MDaN) - 2NxabM'' - LnM„) = 

= {r'^Kf\t,^)+o'^{r'^))Dax\ (2.4.3) 

for some P,Q < —k and if P and Q are equal to —k then k^^^ and k[''^'' have even parity, respectively. 
Note that (2.4.2) and (2.4.3) can always be solved for M and Ma for any given functions K^^\t, ^) and 
i^i^\t,^)- (2.4.2) and (2.4.3) will be called the asymptotic Killing equations. We call the vector field 
K"^ G An asymptotic Killing vector with respect to the foliation determined by the lapse A'' if its components 
M and M" are solutions of the asymptotic Killing equations. Clearly, the notion of the asymptotic Killing 
vectors is less sensitive to the deformation of E and N than that of the strong asymptotic Killing vectors, 
but it is still not independent of the foliation that N and S define. In particular, as we will see below, 
the asymptotic Killing vectors defined with respect to a foliation for which the lapse N tends to zero are 
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different from those defined with respect to ones for which TV — ^ 1 as r ^ oo. The set of the asymptotic 

Killing vector fields will be denoted by . Obviously, C A^ C An for any allowed time axis and 
An can be injected into A. 

Substituting (2.4.1) into the asymptotic Killing equations (2.4.2) and (2.4.3), we obtain 

A = -2(i?ij/3J -pijSj), (2.4.4) 

=2(Bi/3j -A^j) -2(i?ikp''j -piki?''j), (2.4.5) 

and, if E,F < 0, we also have 



T = -2 (Ti /?' - Ti ) , (2.4.6) 
Ti =2(tA -rSi) -2(rVji -r^Ra), (2.4.7) 

independently of k^^^ and . Thus the coefficients Bi , i?i j , T and Tj in both of the elements of A'^ and 
A^ satisfy (2.4.4-7). To calculate the Lie bracket of two asymptotic Kilhng vectors K^, A'" e A^ exphcitly, 
we should compute only the last three terms on the right hand side of (2.1.15), because, by (2.4.2) and 
(2.4.3), the first two are at most of order 0(r^~'^) in general, and zero if K"', K°- e A^^. The last three terms 
are 



+ 4ri-'^^(i?kiBj -i?ki5j)gW-J +o°-(ri-^), 



(2.4.8) 



-(iJkj^m' -^kjiJm')?^*'^-''" -(iikj^-in - ^kj n ) 9^'=^" ' ) + 0°° (r^''^) , (2.4.9) 

(mD^M - MD°-M^DaX^ = 4x'' (Sk^' - SkB') + 2{tB^ - TB') + 

+ 2r«(A(^)B' -^Bk(GA(^)^ + (ajM(^))(og^' -^)))+o°°(r^)- 
-2r^(/.(^)B■ -^Bk(V"^)^+(ajM(^))(o?i' -^)))+o~(r^), (2.4.10) 



where e.g. djfj,^^^ denotes the partial derivative of ^^^\t,^) with respect to its argument fjv'^ = 
^(fc)ij jj^ik ogjigW^ j^j^j ^ggj ^YiaX E,F,G,H < (1 - fc). The right hand side of (2.4.8-10) have the 
form of (the components of) an asymptotic Killing vector, and, as a simple calculation shows, the coefficients 
in K"- := [A", A"]", given explicitly by 



Si =2(i?ij5j -.RijSJ), (2.4.11) 

^ij =2(i?ik^''j -Ri^.B^j -ASj), (2.4.12) 

fi = 2(rj i - TJ iij i + TBi -TBi), (2.4.11) 

f = 2{TiB' -fiB'), (2.4.14) 
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also satisfy (2.4.4-7). However, in general neither nor A'^ form a Lie algebra with respect to the spacetime 
Lie bracket. In fact, for any two spacetime vector fields K"^ and K"^ one has L[j^ K\gab = i-'KX'K9ab ^i^iii^KfJab, 
and for K°-,K'^ G A^ one can form the parts t"t^L[K;,K]5a6, Pci^^[K,K.]gab and f'"-fd^[K,K]5afe- Then, using 
(2.1.1) and the asymptotic Killing equations, one can show that these parts contain terms like MN~'^ or 
MN~^ times factors of order 0{r~''). Thus for general lapse function A'' the order of these parts is not 
0(r~'^). Similarly, A^ does not close to a Lie algebra either. On the other hand, by (2.4.8-10) both A^ and 
A^ are 'essentially' Lie algebras. Next we clarify in what sense do they form Lie algebras. 
Let Q denote the set of the special elements (i/, v"") of A, where 

u{t, x'^ ) = r^!/(^) ^) + 0°° (r^) , lyi {t, x'^ ) = r^i^i^"^^ ^) + (r^) , (2.4.15) 

for some M,N < {1 — k) and the leading terms have odd parity if M = 1 — A; and N = 1 — k, respectively. 
Repeating the construction above, Q can also be considered as the space of the spacetime vector fields 
vt"' + v°', where v is considered to be the lapse of a (maybe degenerate) foliation and t"" is the unit normal to 
the leaves of this foliation. If a lapse N is given on S then we can define Qm in a quite analogous way as we 
did above, and introduce Gf := Af n Gn- Then by (2.4.8-10) behaves as an ideal in Af: [K, kf G Qf 
for any /c" G and K'^ G A^ . We will see in subsections 3.2 and 3.3 below that, at least in the Hamiltonian 
framework, the theory's gauge transformations are generated by precisely the elements of Q. Hence we call 
the elements of Q gauge generators. Since C A^ is a subspace, one may form the quotient space A^ /Q^ . 
This is spanned by the coefficients Bi , , Ti and T, and by (2.4.11-14) it can be endowed with a natural 
Lie algebra structure. Similarly, although A*^ is not closed with respect to the Lie bracket, it almost closes: 
just by (2.4.4-7) and (2.4.11-14) the Lie bracket [K, K]" of any two K", G deviates from an element of 
A^ only by an element of Qn. To determine the structure of Af /Q^ and of A^, we must evaluate (2.4.4-7) 
and (2.4.11-14). 

Applying (2.4.4-7) to N and Ni themselves too, we obtain that r, n , /3i and pij are constant, i.e. 
apart from the gauge generator contents, the coefficients of an asymptotic Killing vector with respect to the 
differential topological background defined by itself are time independent. In particular, r = ri = /?i = Pij = 
corresponds to time axes that are pure gauge generators = vt" + v"', whenever the components T, Ti , Bi 
and of G Af arc all time independent, and the corresponding M and Mi reduce to those given by 
Beig and O Murchadha. However, we saw at the end of subsection 2.2 that such a time axis does not provide 
an appropriate framework in which even the familiar Killing vectors of the Minkowski spacetime could be 
recovered. To be able to recover them we had to assume that, with the notations of the present subsection, 
T = 1. Thus, if r = 1 and Ti = /3i = pij = then the coefficients T, Bi and i?ij are time independent, but 
Ti [t) = Ti — 2tBi , where Ti is constant. Thus, the corresponding asymptotic Killing vector K"^ has exactly 
the same structure as that of the general Killing vector of the Minkowski spacetime in the coordinate system 
based on an asymptotically flat hypersurface approaching the Cartesian one. Therefore, the notion of the 
asymptotic Killing vectors does depend on the fohation coming from S and N. Naturally, by (2.3.5) QmlK'^] 
is still not conserved for general asymptotic Killing vectors, but, as we will see in subsection 4.1, the total 
energy-momentum and angular momentum of the matter -|-gravity system are already conserved. 

If F, ... > then the translation generators T and Ti in (2.4.1) cannot be isolated from the diverging 
gauge generators, and hence by (2.4.11-12) the quotient space Af /Q^ endowed with the Lie product (2.4.11- 
14) is isomorphic to the Lorentz Lie algebra so{l,n). If, however, E,F,... < then T and Ti are well 
defined in (2.4.1). In particular, if E, F, ... < 0, then asymptotically T and Ti dominate the (asymptotically 
vanishing) gauge generators. For E,F,... = the gauge generators do not vanish asymptotically (which 
therefore can be interpreted as supertranslations), but in the limiting case E,F, ... = (1 — A;) = the gauge 
generators have odd parity, while T and Ti , being independent of the spatial coordinates, have even parity. 
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Thus the odd parity supcrtranslations arc proper gauge generators, while the even parity ones, which are 
called the translations, belong to — . Thus, by (2.4.11-14), /Gf is isomorphic to the Poincarc Lie 
algebra. Therefore, the algebra of the asymptotic Killing vectors modulo gauge generators depends on the 
fall-off of the metric: for slow (Q < k < 1) fall-off we have only the Lorentz Lie algebra, but for faster (k > 1) 
fall-off translations emerge naturally and we have a Poincare structure for /Q^ . This result should be 
intuitively obvious: if the asymptotic end is asymptotically flat in any sense then it is becoming spherical 
asymptotically and hence the rotation group (and its relativistic extension, the Lorentz group) emerges 
naturally, but the displacements of the centre of the asymptotic rotations become asymptotic symmetries, 
i.e. the asymptotic translations emerge and hence the symmetry group is the Poincare group, only if the 
geometry falls-off rapidly enough. Finally, we note that the fact that the coefficients Bi , , Tj and T 
satisfy (2.4.4-7) independently of whether belongs to A^ or A'^ implies that the factor spaces A^ /Gf 
and AyG^ are isomorphic, where t/| := Gn n AP^, and hence AP^/G'^ also has the Lie algebra structure that 
Af/G^ does. 



3. The Hamiltonian phase space of the vacuum GR 

3.1 The (partially reduced) phase space and the constraints 

Based on the m = n -\- 1 decomposition, by the a priori configuration variables we would have to mean the 
fields TV, N°' and Qab on a connected n dimensional manifold S of subsection 2.1 (see also [27]). As is well 
known, however, by carrying out the full Hamiltonian analysis with these variables, the fields N and N"" would 
turn out to be pure gauge variables, and the a priori Hamiltonian phase space could be partially reduced to 
the cotangent bimdle T* Q of the (partially reduced) configuration space Q := { qat + boundary conditions }. 
Thus our analysis will be based on the partially reduced configuration space Q and its cotangent bundle 
T*Q. Let qab{u), u € (— e,e), be any smooth 1-parameter family of metrics on S from Q for some e > 0. 
Then we define Sqab '■= {<iqab{u)/du)u=o, which is the tangent vector of the curve qab{u) at qab '■= ?a6(0), i.e. 
Sqab S Tg„i,Q- Obviously, Sqat satisfies the same boundary conditions that qab does. The elements of T*Q 
are the pairs {qab^P""^), where 

P"*" ■Tq^.Q^^- Sqab ^ (r^ Sqab ) := / P^^abd^x. (3.1.1) 

Thus the canonical momentum p""^ is a contravariant symmetric tensor density of weight one on E, which 
is a 1-form on Q, and the requirement of the finiteness of its action on the tangent vectors Sqab gives 
boundary conditions for p"*". The symplectic 2-form on T*Q is the canonical one: for any two tangent 
vectors {SqabtSp"^), {S'qab^S'p"^) £ T^^^^^ pai.)(T*Q) the value of the symplectic 2-form O^^^j^ ^o!,) is 

■.[{5qab,Sp''''),{S'qab,6'r'')) ^ jjySr'S'qab - <5'p"^<5gab)d"a;. ^^'^'^^ 
In terms of the metric and extrinsic curvature the canonical momentum is well known to be 

r' = ;^V^(x"''-X9"'')- (3.1.3) 
Thus the extrinsic curvature (i.e. by (2.1.2) the velocity Qab) can be expressed by the momenta: -\/|^X^^ = 

2l^[P - -(J^ZTTjP ■' Qefq )■ 

The analysis of the field equations in the previous section lead us to the link Z = fc + 1 between the r~'^ 
and r~' a priori fall-off of the metric and extrinsic curvature. Thus, via (3.1.3), we obtain the fall-off 
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P^Hx^) = ;j^P'-'^*^'H^)+o(r-('=+i^) (3.1.4) 



and the parity condition 



pij = (3.1.5) 

r r 

i.e. p'j (*^+^^ is of odd parity. In addition, the 'rest' tt"*" := p"*" — also satisfies the conditions 

o£>c7r"^ = olr-'^-^), o^'do^'cTr"'' = o{r-''-^), ... (3.1.6) 

These imply that oDe,-oDerP"^ = ©(r-^^+^+i)), s = 1, 2, and the parity of the leading term is (-)^+^ 
Following the notations of subsection 2.3, this property of the rest will be denoted by iTab = o°°(r~(*+^)). 
By the requirement of the finiteness of (3.1.1) p""^ must satisfy 

£ p"^5qabdSr = o(r-i) (3.1.7) 

for any Sqab, which implies that n < k + I = 2fc + 1, i.e. k > ^(n — 1). (If wc wanted (3.1.7) to be 
satisfied without global integral conditions, and in particular the parity condition, then we would have to 
require n < k + I, which turns out to be too strong.) Thus the canonical momentum can be interpreted 
geometrically (as the integral kernel of a 1-form on Q) precisely when A; > |(n — 1). In particular, if 
m = n + 1 > A, then fc > 1, and hence the Hamiltonian framework already excludes the possibility of a 
slower, e.g. r~^, fall-off, and yields the Poincare structure for /Q^. Slow fall-off is allowed only in 3 
spacetime dimension. Since Sp'^'' satisfies the same boundary conditions that does, these asymptotic 
and parity conditions ensure that the canonical symplectic 2-form Cl given pointwise by (3.1.2) is also well 
defined. 

As is well known, although the two vacuum constraints c — and Ca — 0, given by (2.1.5) and (2.1.6) in 
terms of the Lagrangian variables qab and (jab, do depend on the lapse N and the shift A''", their expressions 
by the canonical variables Qab and p"*". 



C := c{N, A^^ q.f, q,f{N, (?cd,F'')) = 
Ca := da{N,N',q,f;q,f{N,N',q,a,p''')) = -2qabD,f^, (3.1.86) 



arc independent of N and N°- . Here c := y\q\c and Ca ■= ylql^a, the density weighted Lagrangian con- 
straints, and we used the definition D^tI^- - :— \q\^ DeTI^'" of the covariant derivative of the tensor density 
T^ -; := \q\^Tl^--- of weight w, where Tj^ -- is a tensor field. Asymptotically 

C(?e/,rO = 0(r-('=+2)), Ca{qef,p'n = 0(r-(*=+^)), (3.1.9a,6) 

and the leading terms are of even parity. The constraint functions C : T*Q — > R defining the constraint 
'surface' F in T*Q by their vanishing are 

C[i/,z.°] := J^{c{q,d,p"')'y + Ca{qcd,p"')i^'')d"x. (3.1.10) 
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Here v and //° are smearing (test) fields on S. The requirement of the finiteness of C\v, i/"] yields that the 
smearing fields v and Pa ='■ t^i oDaX^ must satisfy 



z/(t,^'^) =r^z/W(t,=^)+o°°(r^), (3.1.11a) 
z.; (t, x'^ ) = r^i/i {t, — ) + o~ (r^) , (3.1.116) 



where M, A'' < k + 2 — n= {1 — k) + {1 + 2k — n) and if the equality holds in these inequalities then u^^^ and 
respectively, must be odd parity functions. Note that, in particular, for the slowest possible fall-off of 
the metric, k = ^{n — 1), the powers M and N would not be greater than (1 — k). Hence the corresponding 
spacetime vector field fc" = vf^ + v"- would belong to Q of subsection 2.4. 

3.2 The constraint algebra 

Let {qab{u),p°'^{u)) be a curve in T*Q through the point {qab,P'^^) with tangent {6qab,Sp'^'')- Then the 



derivative of the constraint function C[z/, u^] in this direction is 



- ^ D,(^v{q'^^D-5q,^ - q^'^D^q,^) + ^{q'^'D^u - q'^^D-u)5qab+ (3-2-1) 

+ i2v-p''5qab - ^'P^'^hab + 2uJp-^)]^\A^x, 



where 



5C[v, z/' 



Sqab 2k L V \q\ ^ (n - 1) V 

+ D"D''v - q"^'DeD^v^ - ^'^Cq"'' - L^p"^ (3.2.2) 

=^n^'^^ - (^p + ^'^'i'^"- 

Here we used the definition Lxr^' ■ := \q\^ {L^Tl^'" +w{div X)Tj"' ■) of the Lie derivative of the tensor density 
'''b.'.'.' '■— \q\^'^b.'.'.' of weight w along the vector field X'^, where divX is the divergence D^X'^ of X" with 
respect to the natural (metric) volume form. Since we would like to recover e.g. the familiar field equations 
in the Hamiltonian framework in their standard form instead of some of their distributional generalizations, 
we must require the functional differentiability of various functions on the phase space in the strong sense of 
[21]. Thus, in particular, the boundary terms in (3.2.1) must yield zero. Evaluating the leading order and 
parity of the terms in the total divergence of (3.2.1), it is easy to check that the fall-off and parity conditions 
imposed on u and Ui in (3.1.11) already imply the vanishing of the integral of the total divergence. Thus 
Civ, !/°] are already functionally differentiable. 

The analysis of Bcig and O Murchadha given in Appendix A of their paper [5] shows that the vanishing of 
the functional derivatives of C'liy, z/°] together with the constraints C[z/, i'"] = themselves imply the vanishing 
of ly and i^". Therefore, the constraint 'surface' T C T*Q that Cliy,!^""] — defines is 'nondegenerate'. It 
might be interesting to note that in the closed case (i.e. if S is compact with no boundary) Clu, z^"] does 
have critical points on the 'surface' C[y,v'^] = 0, and these critical points represent flat spacetimes [23]. 
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The Hamiltonian vector field of a (functionally differentiable) function F : T*Q ^ R is defined to be 
the vector field Xp on T*Q given explicitly by Xp = {6F/6qab, —SF/6p°'''), and the Poisson bracket of two 
differentiable functions, F and G, is defined by {F,G} := 2Q{Xp,Xg) = Xp{G). Let v, v and f i , Vi have 
the structure (3.1.11a) and (3.1.11b), respectively. Then the 'components' of the Hamiltonian vector field of 
the constraint function C[u,v°'] are given by (3.2.2) and (3.2.3), and the Poisson bracket of the constraint 
functions C[0, and C[0, C*"] is 

\i"x = 

(3.2.4) 



In general, the integral of the total divergence on the right is not zero. The condition of its vanishing is 
N + N < k + 3 — n, and this also ensures the existence of the second integral and that [v, P]"' has the structure 
(3.1.11b). The Poisson bracket of C[0, i/"*] and C[u, 0] is 



{C[0,u-],C[9,0]} = ^J^D,(u{R-f - iRS'^y + {DfU-){Dfu) - {D'^D){Dfvf) 



+ 



(3.2.5) 



The vanishing of the integral of the total divergence can be ensured by N + M < 3 — n. Furthermore, the 

condition of the finitcncss of the second integral is N + M < k + 3 — n, which, at the same time, ensures 
that i'°'Dai> has the structure (3.1.11a). Finally, the Poisson bracket of C[i^, 0] and C[i?,0] is 



{c[z/,0],C[i7,0]} = 2 j Da{up''''DbU-vp''''Dbv)<rx- j Ca{vD''u - vD°'v)A''x. (3.2.6) 

Here the integral of the total divergence is vanishing if M+M < k + 3 — n, which condition, at the same time, 
ensures the finiteness of the second integral and that uD'^v — uD"^P has the structure (3.1.11b). Thus, to 
summarize, in addition to (3.1.11), the orders of the smearing fields must also satisfy N+N, M+M < k+3—n 
and N + M < 3 — n. These conditions can be satisfied by requiring that the powers M and A'^ in (3.1.11) 

satisfy 

M,N<{l-k). (3.2.7) 

By > ^(n — 1) we have 2 + k — n = {1 — k) + 2k + 1 — n > (1 — fc), i.e. the condition (3.2.7) is stronger 

than M. N < 2 + k - n, and, in fact, N + M. N + N,M + M < 2{1 - k) < (3 - n) < 3 + k - n. Note 
that, without further restrictions on the power k or the dimension n, this is the greatest possible bound for 
M and N, because, for the allowed smallest value of the rate of the fall-off of the metric, k = ^(n — 1), 
N,M <{l-k) = 2 + A;-n and N + M <2{l-k) = 3-n. Thus, under the condition (3.2.7), the expressions 
(3.2.4-6) give the familiar Lie algebra C of the constraint functions with the Lie product 

\c\v, z/"], C[P, P"]} = -C[L^i? - Lpi^, [v, - {vD^v - vD^v)]. (3.2.8) 

In particular, the Hamiltonian vector fields of the constraint fimctions are tangent to F on F, i.e. C is a 
so-called first class constrained system. Next we clarify how this constraint algebra is related to Q. 

The fall-off conditions (3.2.7) show that the smearing fields u and v"' are precisely the elements of Q. 
Thus, with the notation fc" := vt"' + u"" in the spacetime picture, we can write CfA;"] := C[u, y""], and then 
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C : G ^ C : k°' 1-^ C[k°'] is surjective, which is obviously Unear. If C[k°-] were the zero function in C for some 
k"^ S G, then the corresponding Hamihonian vector field X(7[fca] would also be zero. Thus, by the result of 
Beig and O Murchadha above, fc" = would have to be held, i.e. C is a vector space isomorphism. But 
in general this is not a Lie algebra isomorphism. If, however, we fix a lapse A'' and fc",^" S Gn, then by 
(2.1.15) and the product law (3.2.8) we have 

{c[fc"],C[fc'^]} = -c[[A:,fc]"] +c[(rt'' + 2g'^'')(i/V(6fcc) -i^V((,fce))i']- (3-2.9) 

In general the second term on the right is non-zero even for asymptotic Killing vectors fc°, fc" G G^ for some 
On the other hand, if we restrict the vector fields fc" and k"" further to be in G^ too, then the second 
term in (3.2.9) is vanishing. Thus although G^ is not closed respect to the Lie bracket, the restriction of the 
constraint function C to G^ mimics the injective Lie algebra (anti-)homomorphisms: the Poisson bracket of 

C[k°-] and C[k°-] for fc", fc" S G^ is just (minus) the constraint function at [k, fc]" € Gn- 

To discuss the linear isomorphism C : G ^ C further, recall that the flow on the phase space generated 
by the Hamiltonian vector field Xc[fe=] = {SC[k'^]/6p°''', —dC[k'^]/dqab) is the congruence of the integral curves 
of the difii'erential equations 



dqab{u) (5C[fc«] dp''''{u) dC[k 



du dp"'" du Sqa 



b 



(3.2.10) 



By (3.2.2), (3.2.3) and (3.1.3), on the constraint surface these are precisely the vacuum evolution equations 
(2.1.2) and (2.1.7) with lapse v and shift with respect to the coordinate time t = u. Let (pu be the 
local 1-parametcr family of diffeomorphisms (on an open neighbourhood of S in the spacetime) generated 
by fc° S G- Then, for small enough u, its action on the canonical variables as fields on this neighbourhood is 
Qab Qab — uQab and p''^ i-^ p°'^ — up""^ . Thus the canonical variables qab and p"*" are changing along the flow 
generated by the Hamiltonian vector field of C[fc"] on F exactly in the same way as under the action of the 
diffeomorphism generated by —k°' on the spacetime. Thus one may say that the flow of Xc[k'^] on F is the 
natural lift of the flow of — fc" from the spacetime to the constraint surface. In the next subsection we show 
that the theory's gauge transformations on the constraint surface are generated precisely by the elements of 

G. 

3.3 The gauge transformations 

Since the Hamiltonian vector fields Xj^j^oj, A;° G G, are tangent to F on F, they belong to the kernel 
distribution 



ker 0|r := {(<5to, 5^') e TF | ^{{dqab, Sp'^'), {S'q^b, S'p^^")) = V {5'qab, S'p'^") € TF } (3.3.1) 

of the pull back to F of the canonical symplectic 2-form O. Since O is closed, this kernel distribution is always 

integrable. But, by definition, the reduced phase space, representing the physical degrees of freedoms, is 
the pair (F, f2), where F is the set of the integral submanifolds of kerf2|r and is the (necessarily well 
defined) projection of O to F. Thus any integral submanifold of ker 0|r through a given point {qab^P"''') € T 
is projected to a single point of F, and hence the vector fields on F belonging to the kernel distribution 
should be interpreted as infinitesimal gauge motions, i.e. generators of gauge transformations. Therefore, 
the Hamiltonian vector fields Xc[k'^] generate gauge transformations on the constraint surface for any fc" G G- 
In this subsection we show that the converse of this statement is also true, namely that any vector field on 
F belonging to kerf2|r (and represented by smooth fields on S), i.e. any infinitesimal gauge motion, is 
necessarily a Hamiltonian vector field Xc[k<^] for some fc" G G- Thus first let us discuss this kernel. 
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By the definition of fl the kernel of Q\r at {qab,P"^) consists of all the vectors (5gfo6, tangent 
to r for which Jj.{Sp''^5'qab - S'p''''6qab)d'^x = for any vector {S'qab,S'p''^) tangent to T at {qab.p"^)- 
To evaluate this condition for {Sqab,Sp°'''), wc should take into account that 6'qab and S'p"''' are not inde- 
pendent. In fact, since T = {(go6,p"'')| C{qef,p''^) = 0, Caiqef,?"^) = 0}, the tangents {S'qab,S'p°-^) of 
r at {qab^P"'^) must satisfy the 'linearized constraint equations' 5'C := {■^C{qef{u),p^^ {u)))u=o = and 
^'Ca '■= {■^Cai'lef{u),p''^ {u)))u=o = 0- Multiplying them by an arbitrary function A and spatial vector field 
A", respectively, and adding them together we obtain 



(3.3.2) 



= A5'C + A«5'C„ = ^-^^5'qab + ^-^^5'^' + Dj^^\5'qa,{q'^''D^X - q^^D^X)- 

Oqab op ^. ZK 

- ^^MK1"'D'^'1-b - q^'^D'S'qab) - 2p-<'\'5'qab + X^p^'S'qab - 2A„5'p-}, 



where SC[X, X^]/dqab and SC[X, X'^j/Sp"^ are given formally by (3.2.2) and (3.2.3). Adding the integral of 
(3.3.2) to J^iSp"'' 6'qab - S'p''''Sqab)d'''x = we obtain 

+ ^-^^y^lab - {Sqab - ^-^^yP"' + DaW'^iX, Xn}d-x = 0, (3.3.3) 

where DaW'^{X, X^) denotes the total divergence in (3.3.2). But from (3.3.3) we can read off the vanishing 
of the coefiicients of S'qab and 6'p'^^ (by the Lagrange lemma of the elementary calculus of variations) only 
if the integral of DaW'^{X, X^) is vanishing. Prom (3.3.3) it follows that it is vanishing if A, A° have the 
asymptotic form (3.1.11), and hence, in particular, C[A, A°] exists. Therefore, the generators X, A" of the 
infinitesimal gauge transformations are precisely the smearing fields. 

Since the vanishing of the Hamiltonian vector field Xc^^.aj implies the vanishing of the vector field 
k" itself, the gauge transformations generated by C[k"'] are effective. The condition ensuring that the 
infinitesimal gauge transformations close to a Lie algebra are the stronger fall-off (3.2.7), which are precisely 
the fall-off properties of the elements of Q. 

3.4 The Hamiltonian 

The aim of this subsection is a concise rederivation of the Hamiltonian of Bcig and O Murchadha, to 
determine the exact and most general boundary conditions for M and Af^ for which the Hamiltonian is well 
defined and differentiable, and to see that M and M" may still be arbitrary functions of time. Thus let us 
start with the 'basic Hamiltonian' 

ffo[M,M"] :=C[M,M"] = j (CM-h4M")d"x (3.4.1) 

and determine that total divergence DaZ" for which the 'total Hamiltonian' H[M,M"-] := C[Af, Af"] + 
Jj, DaZ^d^x is well defined even for the pairs {M, M") defining general allowed time axes, or at least 
asymptotic Killing vectors. Here Z° is expected to be a local expression of the canonical variables, the fields 
M and Af" and their spatial derivatives up to some finite order. 

Let the difference of the physical and background connections bo characterized by T'^j^X'' := (Da — 
oDa)X'', which T^^ is given exphcitly by T^^ = ^q'"^{-oDdqab + oDaqbd + oDbqda)- Then the physical 
curvature scalar of qab can be written as ii = q'^''q'"^{oDa oDbqcd - oDa o-Dc^bd) + qcdq^^q"^ ^Ib^tf - ^le'^tf)^ 
and, following Beig and O Murchadha, we write 
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4^2 1 , 

+Mq,M''V^ {Kb^if - K,Ttf) + Vlid"x. 

(3.4.2) 

The integral of the second and third terms on the right is finite if 

M(t, ) = 2x^ Bu (t) + T(t) + r-^M^-^' (t, — ) + o°°(r-^), (3.4.3a) 

r 

where K < 2k + 2 — n, and if the equahty holds here then /x^^^ {t, ^) must be an odd parity function of 
Here we also had to use k> ^{n—1) if Bi ^ 0. On the other hand, without additional restrictions on k and 
n, the first term has finite integral only for those functions (3.4.3a) in which both Bi and T are vanishing. 

Since the first term is not a pure total divergence, we should write this as the s\im of a total divergence and 
terms that already yield finite integral even for M above. Beig and O Murchadha wrote this 'wrong' term 
as 

Mq''\^^oDa{oDkqcd - oD^Qbd) \/\q\ = 
=oDa(Mq''\<^{oDbqod - oDoqbd)V\Q\) ' MqC^^^^^'^v^) (oA^cd - o£'c5fed)- 

- {oDaM)q''^q'"^(oDb{qcd - o9cd) - o-Dc(%d - oqbd)) VH = 
=oDa(^Mq''''q'%oDbqcd - oD^qbd) VWl) - MoDa{q''''q''' VW\) {oDbqcd " O^cgw)- 

- oi?a((oAM)q"V''('Zcd - oqcd)VW\- {oDbM)q''-q''''{q,d - o9cd) /M) + 

+ oDb[{oDaM)q-y^\j {q,i - oqcd) - oDc[{oD,M)q'^'>q'''^) [q^ - oqbd) = 
= - Da(^Mq''Y'^{oD,qM - oD^qcd) + {oD^M)q''\-\q,a - oqcd) - {oDcM)q''\-'^{qM - oqbd))^M+ 
+ oDa oDbUi^q'^^-^q^d - oqcd) - q''V{qcd - oqa)) V\q\ + 

+ oDa{q''''q'"^V\q\) (-M{oDbqcd - oDcqbd) + (oDbM) {qcd - oqcd) - [oDcM) {qM - oqbd)) ■ 

(3.4.4) 

The integral of the second and third terms on the right of (3.4.4) is finite if M has the form (3.4.3a) where 
K satisfies the stronger condition K < k + 2 — n, and if the equality holds in this inequality then ij,^^\t, ^) 
is an odd parity function of Then, since K<k-\-2 — n<2k + 2 — n, 

H[M,Q] := C[M,0] [ D A M q'^\'\oD ,qM - oDbqcd) + (oAM)(?'^^^'^(ged - o<?cd)- 

^ (3.4.5) 

- (o£'cM)g«V''(96d - 096d)}dS 

is already well defined. If {5qab, Sp"'^) is any tangent vector at {qab,P"'^) G T*Q, then the derivative of H[M, 0] 
in the direction {5qab,Sp"'^) is 

5H[M,0] = l^(^^£^5q.b+^-^^Sp'^')d-x-^ i?„{Mg'^^^'^V^(r^,^g,e - r^6'^9ce) + 

+ S{q''\'"^^M) (M{oDeqbd - oDbqcd) + (o-DfeM) (ged - oqcd) - (o^'cM) {qbd - o5fed)) }d"a;, 

(3.4.6) 
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where 5C[M,Q]/5qab and 5C[M,0]/ dp"^ are given formally by (3.2.2) and (3.2.3). (Strictly speaking, these 
are not functional derivatives of C[M, 0], because the constraint function C is not well defined for M above.) 
Since, however, the integral of the total divergence on the right of (3.4.6) is vanishing for the functions M 
given by (3.4.3a), H[M, 0] is functionally differentiable with respect to the canonical variables too. 
We can write 

Ho% M«] = -2 ^ (D„r^)M6d"a; = 2 ^ (p''\D(,Mk) - r'>Tl^,M, - Da{r'>Mb))d^x. (3.4.7) 
The integral of the first two terms on the right is well defined even for vector fields M° of the form 



Mi {t, x*^ ) = 2x'' Rki it) + Ti (t) + r'^ixf^ (i, — ) + o°° (r^), (3.4.35) 



r 

where L <k ^2 ~ and if the equality holds in this inequality then /ij (t, ^) is an odd parity function 
of Note that, to prove the existence of the integral of the first two terms for nonzero i?ij in (3.4.3b) we 
also had to use k > ^(n — 1). The integral of the third term on the right of (3.4.7) is, however, finite only 
for Rki{t) = and Ti (t) = 0. Thus 

i?[0,M"] := C[0,M"] +2 / £)„(p"''Mb)d"a; = [ p'^^hMQabd^x (3.4.8) 

is well defined even for vector fields M° with the structure (3.4.3b). The derivative of H[0,M°'] in the 
direction {SqabtSp"^) is 

6H[0, M'^] = £ (hMQab Sp'^' - Lm^' Sqab + D,{M-p''''5qab))d"x. (3.4.9) 

Since for the vector fields above the integral of the total divergence is zero, H[0,M'^] is differentiable 
with respect to the canonical variables too. 

Therefore, the Hamiltonian i?[M,M°] := H[M,0] + H[0,M"] of Beig and 6 Murchadha is finite and 
functionally differentiable with respect to the canonical variables even for M and given by (3.4.3) with 
K,L < fc + 2 — n, and if the equality holds in these inequalities then /j,^^^ and /i^'' must be odd parity 
functions of respectively. However, the spacetime vector field K"^ := Mt"" + M" is still not needed to be 
an asymptotic Killing vector field with respect to some foliation (and not even an allowed time axis), because 
the fall-off rates K and L are still required only to satisfy K,L < k + 2 — n (instead of K,L < {1 — k)). 
Moreover, i?i j , Si , Ti and T may still have arbitrary time dependence. 

3.5 The algebra of the Hamiltonians and the asymptotic symmetries 

On r the system of equations 

dqab _ 6H[M, M^] d^_ _SH[M^M^ 

dt Sr' ' dt Sqab ' ^ ^ ^ 

defining the Hamiltonian flow on F, is precisely the system of the vacuum evolution equations with lapse M 
and shift M". However, this system still does not preserve the boundary conditions (2.3.1a-3a) and (3.1.4- 

6) for the canonical variables, because the regularity and functional differentiability of the Hamiltonian 
H[M, M"-] implied only K, L < k + 2 — n, which is weaker than K,L < (1 — k). Thus, based on the analysis 
of subsection 2.3, we must require that the powers K and L satisfy 

K,L<{l-k). (3.5.2) 
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Therefore, = Mt°- + A/" already has the asymptotic form (2.4.1), i.e. K'^ must be an element of A, and 
the restriction of the Hamiltonian H to the pure gauge generators fj, and fx"' (i.e. for which fxt"" +/Lt'* e S), is 
just C[/x,/x"]. 

Repeating the analysis of subsection 3.2, one can show that (3.5.2) ensures the vanishing of the boundary 
terms appearing in the calculation of the Poisson brackets of two Hamiltonians, H[M, M"] and HIM, M"], 

and, by (2.4.8-10), that M^DaAf-M'^D^M has the structure of a lapse and [M,M]'^ and MD^M-MD^M 
have the structure of a shift satisfying (3.5.2). The resulting Poisson algebra of the Hamiltonians [5] (see 
also [4,28]) is 



[M, M"] , H [M, M"] ] = -H LmM - LmM, [M, M] " - {MD^M - MD^M) . (3.5.3) 

For pure gauge generators (3.5.3) reduces to (3.2.8), and the Poisson bracket of iJ[M, M"] and C[i',p°'] is 
also a constraint function. Thiis the Poisson algebra C of the constraint functions is an ideal in the Poisson 
algebra H of the Hamiltonians parameterized by the allowed time axes if G A. The structure of H can 
be determined easily by considering the Hamiltonians parameterized by the special allowed time axes like 
K'' = 2x^Bk{t)t°', K°- = 2x^ Rj' {t){-£rT, ■■■ etc. With this parameterization of the Hamiltonians (3.5.3) 
shows that, for each fixed value t of the coordinate time, the factor of Ti. with the ideal C is just the Poincare 
algebra. Thus H/C is infinite dimensional. If, however, the coefficients Bi , , Ti and T are restricted to 
be the coefficients in the asymptotic Killing vectors with respect to some as in subsection 2.4, then the 
whole factor H/C would be finite dimensional, and, in fact, the Poincare algebra. 

As we noted in subsection 2.4, the space A of the allowed time axes docs not form a Lie algebra with the 
natural Lie bracket in general. Hence the Poisson algebra 7i of all the Hamiltonians, indexed by the elements 
of A, does not seem to be connected in a natural way to some naturally defined Lie algebra of spacetime 
vector fields. However, restricting the spacetime vector fields K" and K'^ to be from the subspace A'^ of the 
space A^ of the asymptotic KilHng vectors for some and writing HlK"] := H[M,M"], by (2.1.15) the 
product law (3.5.3) takes the remarkably simple form 



(3.5.4) 



If if" and if" are allowed to be from A^ , then the first two terms on the right of (2.1.15) give a constraint 
function with uncontrollable generators, as in (3.2.9), and we have only 



{ii[is:''],ii[if"]} + ii[[ii:,if]"] GC{g)=c. 



(3.5.5) 



Therefore, the set := H{A^) of the Hamiltonians parameterized by the asymptotic Killing vectors K" 
with respect to is the Lie product preserving image of A^ modulo constraints, and on the elements of the 
subspace A'^ the Hamiltonian H preserves the Lie product. 

Finally, let M, N, and M" be as in (2.4.1), and calculate the total time derivative of iJ'[M, M°] 
along := Nt"" + N"". Since M and M" may depend on t, the derivative consists of two terms: 



—H[M,M''] =if[M,M*] + 



6H[M, W 



H [M, M'^] + I if [N, N''] , H [M, M"] | 



6H[M, 
5qab 



(3.5.6) 



= H[M + h^N - LnM, + ND^M - MD^N - [N, M]"] . 

Here first we used (3.5.1) (which, on the constraint surface, are the vacuum evolution equations), and then 
(3.5.3). If e A^, then by (2.1.13) and (2.1.14) the right hand side is zero, while for e Af the right 
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hand side is a constraint function by the asymptotic Kilhng equations (2.4.2) and (2.4.3). Therefore, the 
Hamiltonian of Beig and O Murchadha is constant along any allowed time axis ^'^ modulo constraints for the 
asymptotic Killing vectors if" G , and it is strictly constant for vectors K"^ that are strong asymptotic 
Killing with respect to the time axis 

To summarize: First, to ensure that e.g. the symplectic 2-form be well defined, in addition to the result 
Z = fc + 1 of the analysis of subsection 2.3, wc had to assume that k > ^(n — 1). Then the constraint functions 
are well defined, functionally diffcrcntiablc and close to a Lie algebra precisely for those smearing fields that 
correspond to the elements of G itself. The constraint function preserves the Lie product of the elements 
of the space t/| in the Poisson algebra. The Hamiltonian of Beig and O Murchadha, parameterized by the 
allowed time axes K'^, are finite, functionally differentiable, close to an infinite dimensional Lie algebra, 
and the Hamilton equations preserve the boundary conditions for the canonical variables. The Hamiltonian 
preserves the Lie product of the elements of the space A'^, and preserves the Lie product of the elements 
of the space A^ modulo constraints. The Beig-0 Murchadha Hamiltonian H[K"-] is constant in time with 
respect to ^" for any if" G A'^, but it is only constant modulo constraints for the asymptotic Killing vectors 
if" e Af. 

4. The ADM conserved quantities of matter+gravity systems 
4.1 The ADM conserved quantities 

In the complete Hamiltonian description the matter fields would have to be included. However, a detailed 
Hamiltonian analysis of the matter fields is not needed if we are interested only in the ADM energy- 
momentum and angular momentum, because the value of the Hamiltonian of the matter fields on the 
matter constraint surface (if there is any, as in the spacial case of Yang-Mills fields) is expected to be 
Qm[K"'] for some allowed time axis if" G A. Then the gravitational constraint, i.e. a part of Ein- 
stein's equation, is Qm[K°'] + C[if"] = 0. Thus, the value of the total Hamiltonian on the constraint 
surface, i?[if"]|r + QmlK""], is given by the same surface term as in the vacuum case, and it has the 
structure iJ[if"]|r + Qm[K''] = T{t)p° + Ti{t)p' + i?ij (t)j'j + 2Bi{t)j'°, where T{t), Ti{t), i?ij (i) and 
Bi (t) are the functions appearing e.g. in the form (2.4.1) of the allowed time axis if". If, however, 
if" G Af for some ^" G A, whenever the functions T{t), Ti{t), Rij{t) and Bi{t) can be represented 
by the ^m{m + 1) independent parameters T, T; , i?ij and i?; , the evaluation of ii[if"] + Q,n[K'^] on 
r above defines a linear mapping Af /Gf « A^/G'^ R, whose components P°, P' , J'j and J'° de- 
fine the energy, the linear momentum, the spatial angular momentum and centre-of-mass, respectively, via 
iJ"[if"]|r + QmiK"] =: TpO -h TjP' -|- iijj J'J -|- 25; J'". In particular, if if" G Af and ^" G G, whenever the 
functions T{t), T; {t), Rij {t) and B^ (t) are all constant, then, as we will see, P'^ and P' are just the familiar 
ADM energy and linear momentum, respectively, J'^ is the angular momentum of Regge and Teitclboim and 
J'" is the centre-of-mass given by Beig and O Murchadha. If, however, if" G Af but ^" — t" G G, whenever 
T{t), iiij {t) and Bi {t) are still constant but Ti (t) changes in time as Ti [t) =Ti — 2tBi (just according to the 
expression of the n + 1 form of the boost Killing vectors of the Minkowski spacetime), then P*^, P' and J'-" 
remains the same but the centre-of-mass J' " deviates from that of Beig and O Murchadha by the term iP' . 
Since P- and J-- , a , 6 = 0, 1, n, are elements of the dual space of Af /Gf , it is easy to check that under 
a Lorentz transformation x- x-Ab- of the Cartesian coordinates x- — {t,x^) the energy-momentum P- 
transforms as a Lorentz vector and the angular momentum J-- as an anti-symmetric tensor: P- i— » P-A;,- 
and J-- I— > J--Ac-Ad-, while under a translation x- i-^ x- +ri^ of the Cartesian coordinates P- remains 
intact and J-- i— > J^^ -|- 2r]^P-h Thus under a Poincare transformation of the Cartesian coordinates P- 
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transforms like the energy-momentum vector and J-- as the angular momentum tensor of a Poincare in- 
variant system. We emphasize that to derive these transformation properties the centre-of-mass expression 
of Beig and O Murchadha had to be completed by the term iP' . It might be worth noting that in the 
so-called field formulation of general relativity on a given background [29,30] the centre-of-mass expression 
also contains the extra term tP' . 

The general expression of these quantities in terms of the metric and the extrinsic curvature for any 
allowed time axis K°' is 



Q[i^«] := + = ^ Da{Mq''\^'^{oD,qM - oDtqcd) - {oDbM)q''\^\q,d - o9cd) + 

+ {aD,M)q'^\'''{qM - oqbd) - 2Mt,{x'"' - X(?'")}dS. 

(4.1.1) 

Therefore, we can define the energy-momentum and angular momentum of any asymptotic end by the 
surface term of (4.1.1) even in the presence of matter fields, independently of any symplectic or Hamiltonian 
structure or phase space. The only requirement is its existence, and we assume only that the boundary 
conditions obtained from the investigations of the evolution equations in subsection 2.3 hold. Apparently, 
for asymptotic boosts (i.e. for Bi ^ 0) and rotations (i?ij ^ 0) (4.1.1) gives finite value only if fc = n — 1, 
and for asymptotic translations (i.e. for B\ = 0, = but T 7^ or Ti 7^ 0) only if A: = n — 2. However, 
it is well known that the ADM energy-momentum is finite and well defined even if the metric falls off only 
slightly faster than r~5("-2)^ because the contribution of the slow fall-off 'part' of the metric and extrinsic 
curvature to the ADM energy- momentum can always be written as a constraint: the whole Hamiltonian 
expressed as a volume integral is finite even for the slower fall-off asymptotic ends [10]. We show that, by 
the same reason, the angular momentum and centre-of-mass can be finite even for metrics with r^i("-^^) 
fall-off. In particular, in m = 3 -|- 1 spacctime dimensions the a priori 1/r fall-off of Beig and O Murchadha 
is the weakest possible for which, in general, we can have finite angular momentum. 

To determine the weakest possible power- type boundary conditions coming from the finiteness of Q[K'^], 
let us rewrite it as an integral on S by the very definition (3.4.5) and (3.4.8) of the Hamiltonian: 

Q[K'^] = f {fMa + MM)dS + i / (x"" - X?"') {oD^aMb) - K,M,) dS- 

-^jjyoDa oDkM((f\''\q,a - mod) - q^U'^'iqcd - o9cd)) 

+ i^Da{f\'''^\) [{oDbM) {q,d - oQcd) - {oDcM) {qM - oqbd)) - 
- MoDa{q''\'''^M) {oDbqcd - oDcqhd) + 

+ Mq^M^'q^^'iri.rif. r:,Ttf)VW\ + M{x' - xabx'"')VW\}d''x, 

where M and Ma have the form (2.4.1) for some (unspecified) powers E and F. Suppose that the energy 
density and the momentum density of the matter fields satisfy the fall-off and parity conditions of subsection 
2.3. If Bi 7^ then the condition of the existence of the integrals involving M is < (1 — A;) and A; > i(n — 1), 
and if the equality = (1 — A;) holds then iJ-^^\t, ^) of (2.4.1) has odd parity. Thus, in particular. A; > 1 
and E < must hold if n > 3. If Si =0 then the rate k of the fall-off can be reduced. In fact, the condition 
of the existence of the integrals involving M is E < —k and k > i(n — 2), which, for n = 3, gives the 
well known results k > ^ and E < —k < — i. Similarly, if Rjj ^ then the condition of the finiteness of 
the integrals involving is A; > \{n — 1) and -F < (1 — A;), and if the equality F = (1 — A;) holds then 
p^P(t, ^) of (2.4.1) has odd parity. If iiij = then the fall-off may be slower: A; > \{n - 2) and F < -k. 



(4.1.2) 
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Therefore, the eneryy-m,om,entum, and (relativistic) angular momentum are finite for general time axes K"' 
precisely when k > i(n — 1), but for the slower fall-off k > i(n — 2) the finiteness of the energy-momentum 
is not guaranteed by Bi =0 and = alone, E, F < — fc must also be required. 

This motivates us to consider, for some q < {1 — k), the special time axes K"^ = Mt"" + M° e A with 
the asymptotic structure 

M = T(t) + r Vn*' — ) + 

(4.1.3) 

Mi ^ Ti (t) + r^/ii (^) (t, ^) + 0°° (r^) , i?, F < g. 

A simple calculation shows that they do not form a Lie algebra. If, however, they arc assumed to be 
asymptotic Killing vectors too (whenever the powers P and Q in the asymptotic Killing equations (2.4.2) 
and (2.4.3) should be required to satisfy P,Q < q — I, and T{ and T are necessarily constant), then for 
certain values of q they form a subspacc in which behaves like an ideal of a Lie algebra. In fact, if qT^ 
is the set of the asymptotic Killing vectors = Mf" + Af° whose components satisfy (4.1.3), then by a 
calculation similar to (2.4.8-10) shows that the Lie bracket of iC" G qT^ and K"^ G A^ (where the latter is 
given by (2.4.1)) contains terms of order r~*'. Thus the Lie bracket operation preserves the index q of 
and the components of [K,KY have the structure (2.1.3) provided q > —k. The quotient /qT^ n 
is isomorphic to R"^ and inherits a commutative Lie algebra structure from A^ /Q^ . Thus qT^ may be 
interpreted as the space of the 'g fast fall-off' asymptotic translations in A^ where — fc < g < (1 — fc). 
These asymptotic translations can be singled out even if A; G (0, 1), whenever A^ /G^ is only the Lorcntz 
Lie algebra rather than the Poincare one. The results of the previous paragraph show that for fc > ^(n — 1) 
the space of the translations could be any of gT^, q e [— fc, 1 — fc], but for fc > ^(n — 2) it is just the space 
-kT^ whose elements yield finite energy-momentum. 

For the sake of logical completeness one should note that, strictly speaking, the standard expression 
for the ADM energy-momentum and angular momentum (including the Beig-0 Murchadha centre-of-mass) 
differs from that given by (4.1.1). However, it is easy to see that (4.1.1) coincides with the standard one. In 
fact, for example the first term of the integral (4.1.1) can be written as 

^ Da[Mq'^\'=''{oDoqbd - oA^cd) } ^d"x = 
= j o£>a{M(og"*' - r-^q^^'i''^ + o°°(r-*)) (og^'' - r'^q^^^'"^ -{■ o°°(r-*)) x 
X {oDcqbd - oDbqcd) ^1 + r-'=d + o~ (r-'=) } vToaId"a; = 
= -^lim £ (Moq'^{oDcqad - o£>agc<i) + r--'=F„) ot^V"-M5, 

where (jfC")"^ := o?"^ og'"^gcd\ d is a smooth function, qI'" is the outward directed ogafe-orthogonal unit normal 
to the large sphere Sr of coordinate radius r, AS is the unit sphere volume element, and the 1-form Fa is 
defined by the last equality of the integrands. If K°- e A and fc > |(n — 1), then Fa qv"" = 0{r~'') holds and 
its parity is odd, while if has the structure (4.1.3) with q = —k and fc > |(n— 2), then Fa ov"' = 0(r~^~^). 
However, in both cases the r ^ oo limit of the integral of the term r"~'^~^Fa qv'^ is zero. Similarly, all the 
remaining terms in (4.1.1) can also be written into the form being linear in the physical metric qab and the 
extrinsic curvature Xab, yielding the familiar expression given in [5]. 

Finally, calculate the total time derivative of along any allowed time axis where K'^ e An- 

However, in the present calculations we cannot use (3.5.1), because they are the vacuum evolution equations 
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in the Hamiltonian phase space. In the presence of the matter in the spacetime we must use (2.1.2) and 
(2.1.7). They, the definitions and formulae (3.1.3) and (3.2.2) imply that 



2"" \2k 

Then by (2.3.5), (3.5.3), (4.1.4) and the definitions we have 



XcdX 



2 
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(4.1.4) 



=Qm [M + Lm^V - LnM, + ND'^M - MD^'N - [N, M^] +H[M, M""] + 

+ [h[N, N^] ,H[M, M^]}+J (i- (i? + x' - XabX"') - {Mx + £>cM^)7VdS = 

=H[M + Lm-/V - LnM, M*" + ND^'M - MD^'N - [N, M]"] + 
+ Qm[M + Lm-/V - LisrM, M'' + ND^M - MD^N - [N, Mf] + 

+ ^ (i? + - XafeX"') - m) {Mx + D,M-) TVdE. 
Taking into account the Lagrangian constraints (2.1.5) and (2.1.6) we obtain 

[M, M«] =Q[M + hmN - LnM, M^ + ND'^M - MD^'N - [N, M]«] . (4.1.5) 

By (4.1.1) the right hand side is an (n — l)-sphere integral at infinity with the generators M := M + — 
LnM and M^ := M^ + ND^M - MD'^N - [N,M]''. If K'' is a general asymptotic Killing vector, then 
by (2.4.2) and (2.4.3) the order of these generators is at most 0{r^~''), whenever their parity is odd and 
k > ^{n — 1). If is an asymptotic translation from -kT^ and k > ^(n — 2), then the order of these 
generators is at most 0{r^'^). However, in both cases the right hand side of (4.1.5) is vanishing. Therefore, 
the quantities Q[K°-] are constant in time for the asymptotic Killing vectors K" e . This result is 
analogous to the fact that the Qm{K'^] of subsection 2.3 is constant in time for any Killing vector K"^ of the 
spacetime. Nevertheless, while for the conservation of Qm[K'^] built from the matter field variables K"^ must 
be a genuine Killing vector, the conservation of the analogous quantity of the matter+gravity system 

is ensured even by the asymptotic Killing fields too. 

However, if the time evolution is defined by the allowed time axis but K'^ is an asymptotic Killing 
vector with respect to another i.e. K'^ G , then in general (5[/i'"] is not conserved. In particular, if 
represents pure time translation at infinity (i.e. its lapse part tends to 1), but K"^ G A^ for some G Q 
(whenever Bi{t), (t), Tj (t) and T{t) of K"^ are independent of t), then Q[-ftr''] is not constant in time 
with respect to For example, the time derivative of the centre-of-mass of Beig and O Murchadha with 
respect to above is not zero, that is just the spatial (linear) momentum. 



4.2 The background-independence of (5[^"] 

Let ov"' be the outward directed o9a6-unit normal to the coordinate spheres in E, and let 7 be an integral 
curve of ov"' form some Sr^ to Sr- Then the length of 7 in the physical metric Qab is 
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R= f \l\qab^v^\dr' = j J 1 - ^q'^^ oV' ovb + o{r'-'')dr' = 



r-ro + Aln^+B + o{r-°), if/c=l; 
r-ro + Ar-^+i +B + o(r-'=+i), ii k ^ 1 



(4.2.1) 



for some constants A, A, B and B. This implies, in particular, that 

J_ _ ]_ ^ i o{r-^), if/c = l; 
Rk \0(r-2'=), iffcT^l. 



(4.2.2) 



Therefore, in the definitions (2.3.1)-(2.3.3) of the asymptotic flatness the radial distance r can be substituted 
by the physical radial distance R without changing the structure or the leading terms of qij and Xij • 

To clarify the potential ambiguity both of the notion of asymptotic flatness and the quantities (3[-fl''*] 
coming from the non-uniqucncss of the background metric oQab, let (S, Qat, Xab) b<-> (k, /)-asymptotic;ally flat 
with respect to o^abj and let ofjaft be another backgroimd metric, being flat on E — K . (Without loss of 
generality we may assume that the domain of the flatness of both o(?a6 and oQab coincide.) Thus there exists 
a diffeomorphism (j): 'E — K^T, — K such that o9o6 = 0* o^ofc- For the sake of simplicity suppose that is 
homotopic to the identity Idls-x, i-e. for some one-parameter family ^„ of difl[eomorphisms = Id|s-if 
and = (j). Then we can form the one-parameter family of flat metrics oqab{u) ■= (j)^ oQab on S — if (which 
can obviously be extended to the whole S as, in general curved, negative definite metrics). If y is the vector 
field on S — K generating and its components in the coordinates {x^ } are defined by VoSab ='■ qD^x^ , 
then 

^09a6 := (^09a5(w))|„=0 = Lvo9a6 = (o A ^ + QDjVi^oDaX\DbX^ . (4-2.1) 

Writing Vi in the form 

Vi {x^ ) = 2x'^pki + Ti + r«y/-^^ (^) + (r^'^) (4.2.2) 

for some power R, where the first two terms together is just the kernel of the fiat Killing operator o-D(i ) 
for o^ab) we have 

oZ)(iV-j) =r«-i(i?ot'''4(iV'jf + {d^vl^^){5^)~h)r^ ov"" ov"")) + {r^-'). (4.2.3) 

Its leading term has even parity iff V^^^\^) has odd parity. Since, for sufficiently small u, one has gab — 
oQabiu) = qab - oQab - ioqab{u) - oQab) = Qab - oQab " 5 oQabU + Oiu^), and hcncc, in the o9a6-Cartesian 
coordinate system {x''}, gij - o^ij (w) = Q'ij - oaij -<5oftj u+Oiu^) = r'^q^") +o°°{r^'')-2oDi^iV^)U+0{v?) 
holds. The one-parameter family of coordinate systems {S'(u)}, defined in terms of the coordinates {x^} 
by x^{u) := (l)\^{x^), is Cartesian with respect to the one-parameter family of flat metrics oqabiu), i.e. 
o9mn (m) = — (^mn ■ (To sec this it is enough to recall the definition of the pull-back oQabiu) of the metric 
oQab along (f)u in the coordinates {x'}, viz. o9ij (w) := ^li-^o^mn = -(^mn ^r-f^r^ and to compare 
this with the transformation law o9ij (w) = aq-mniu)^^^^ of the components oqij{u) and o^mn('«) of 
oqab{u) in the coordinates {a;'} and {a;™(u)}, respectively.) Then for sufficiently small u we have x^ = 
x^ + V^{x^)u + O(u^). Thus the components of the physical metric and the extrinsic curvature in the 
o9a6(u)-Cartesian coordinate system {x^{u)}, respectively, are 
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. X . (4.2.4) 

-2«r«-i-'(i?ot;'"<5^(iX?jll^^^^'' +xL'(i^J)^^''^'' - o«'"(a™l^^^^'^)xL'|i^j). o«') + 

Therefore, i/ie asymptotic end {^,qab,Xab), which is {k, I) -asymptotically flat with respect to the background 
metric oQab, remains {k, I) -asymptotically flat with respect to the new flat metrics oQabiu) — (/i^oSah for 
sufficiently small u precisely when R< (1 — fc) and has odd parity for R=l — k. These changes of the 
background metrics wih be called allowed. (Apart from a rigid Euclidean rotation, the non-trivial leading 
terms of gij and Xij &re invariant with respect to the change o^ab '-^ oQabiu) of the background metric iff 
i? < (1 — A;).) Therefore, the corresponding F"'s are special, purely spatial asymptotic Killing vectors. 

The change of the background metric oQab yields a change of the allowed time axes. By (4.2.2) for the 
infinitesimal change of M and Ma given by (2.4.1) one has 



=2x'' (2pkj B-" ) + 2Tj +2r^ot^''pk'(5ii/(^)) +2r^S'y/^^ +o~(r^) +o~(r«), 

= (2a;''(2pkji?ji -2iikVji) + (27^ iiji - 2TVji) + 2r^( ot^'^Pk^ (Sj - t/f V i) + 

+2r«o^^''i?kj (iJo^^^^miV/^^ + (a„y/''^)((5r -ov"^ ot^" ^ni )) +o°°(r^) +o~(r«)) oDax\ 

(4.2.5) 

Thus the allowed change of the background metric acts on K"' = Mt°- + M"- as the diffeomorphism generated 
by the asymptotic Killing vector V: for E,F < [1 — k) the structure of SM and SMa is similar to that of 
M and Ma, respectively, if i? < (1 — A;) or if J? = (1 — k) and v/^^ has odd parity. If, however, Bi = 0, 

= and E,F < -k, i.e. K'' e -kT^ , then the structure of 5M and SMa will be similar to that of M 
and Ma, respectively, only if i? < —k. 

To calculate the change of QlK""] under the allowed change of the background metric let us form 
Qu[K''{u)] by using the one-parameter family of the flat background metrics oqab{u) in (4.1.1) instead of 
o^ab- Then, with the notation SK"- := SMt" + q"'^SMb, a straightforward calculation gives 



= Q[SK-] + — Da[M{Soriq''' -SoK.q"') + {qV - q^'q'") {oDbM)6 oqcdjdi:, 

where 5 qF^^ := i o9°''(~o-Dd^ oQbc+oDcS oQdb+oDbS oqcd)- Thus the spatial momentum and angular momen- 
tum depend on the background metric only through their generator Ma, but the energy and centre-of-mass 
may be ambiguous as a consequence of the non-vanishing of the integral on the right of (4.2.6) too. (Using 
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the expression (4.1.2) for instead of (4.1.1) it can be shown that 5Q[K°-] is finite if i? < (1 — k), or if 

i? = (1 — fc) and V"/^' has odd parity.) First suppose that if" = Mt" + A'/" is an asymptotic KiUing vector, 
and hence E,F < (1 — k). Then by the results of the previous subsection k > ^(n — 1) must hold, and by 
(4.2.4) and (4.2.5) i? < (1 — A:) must be required. Thus (5[(5if''] depends only on the first two terms of 5M 
and SMa in (4.2.5), because the remaining terms are pure gauge generators. However, by (2.4.11-14) this 
is nothing but the transformation law of the energy, and the components of the spatial momentum, spatial 
angular momentum and ccntre-of-mass under the Euclidean transformation coming from the diffeomorphism 
generated by V"". To rule out the ambiguities in the expression of the energy and the centre-of-mass, we 
must require the vanishing of the integral in (4.2.6). Its vanishing can be ensured by requiring R < {2 — n), 
OT R = (2 — n) and that V^^^ be odd parity functions. Next suppose that K'^ e -k^^^ , i.e. = 0, 
B\ = Q and E,F < — fc. Then by the previous subsection Q[K'^] is finite even if fc > ^{n — 2), and by (4.2.5) 
R < — fc, implying that (^[(Jif"] is finite and describes how the components of the spatial momentum trans- 
form under the Euclidean transformation coming from the diffeomorphism generated by V"". (The energy 
remains intact.) The vanishing of the integral in (4.2.6) can be ensured even by < (3 — n). Therefore, 
Qf-ft'"] is unambiguously defined for K'^ e if R < {1 — k), R < {2 — n) and in the case of the equality, 
R = (2 — n), V^^"^ has odd parity; and for K'^ £ -kT^ if R < —k and R < (3 — n). In particular, in m = 3 + 1 
spacctime dimensions and for fc = 1 the angular momentum and centre-of-mass arc well defined provided the 
diffeomorphisms connecting the background metrics tend to the rigid Euclidean transformations like 0(r~^) 
with odd parity generator, or faster. For the condition ensuring well defined energy and spatial momentum 
we recovered the known results of [8,11]. Namely, writing fc in the form fc = i -|- 5 for some d > 0, the energy 
and momentum are well defined if the diffeomorphisms tend to the Euclidean transformation like 0{r~^~^). 

Appendix: Global integral conditions at the null infinity 

Let us consider global quantities associated to the global state of the matter fields in the Minkowski 

spacctime at a given retarded time U, i.e. if the spacelike hypersurface of subsection 2.2 extends to 
the future null infinity Z+ such that its intersection with is the U = constant cut. However, in 
the present context it seems more comfortable to choose this hypersurface to be the null hypersurface 
Afu ■= {(r, X')| T — R = U} instead of a spacelike (e.g. a hyperboloidal) one. Here (r, X'), i = l,...,n, 
are still the Cartesian coordinates introduced in subsection 2.2. For the null Mu the flux integral of T°-^Ki, 
takes the form limi^^oo /g §gT°-^laKbdSR''^~^AR', where la ■= Va(T — i?) = t,, + ?'„, a null normal to JVu- 
Thus if, for the sake of brevity, wc define piU, R, ^) := taT^Ht = /i([/ + R,X^ )+ f{U + iJ, X^)va and 
p' {U, R, ^) := KlT^Hb = K'l {j^iU + R,X^) + a'^^U + R, )vb), then the global energy, spatial momen- 
tum, angular momentum and centre-of-mass measured at the retarded time U at future null infinity are finite 
precisely when +/9d5 = o(i?-"), §g+p'dS = o(i?-"), -p^'X^^j^dS = o(i?-'") and p' p^)dS = 
o{R-^), respectively. These conditions could be satisfied if p{U,R, ^) = ■^p^"'HU, ^) + o(i?-™) and 
p^ iU,R,^) = ^(+ip(U,R,^) + ^-p^"'\U,^)) + o{R-"'), where -p('")(C/,^) is the odd parity part 
of p^"^\U, ^) and ~^(p{U, R, ^) is an arbitrary function with even parity, ^i^p^™) contribute to the energy 
and spatial momentum, but do not to the angular momentum and centre-of-mass. Thus we may call them 
the BS mass aspect of T"*". 
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